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ABSTRACT 

The RCS 2319+00 supercluster is a massive supercluster at z = 0.9 comprising three optically 
selected, spectroscopically confirmed clusters separated by < 3 Mpc on the plane of the sky. This 
supercluster is one of a few known examples of the progenitors of present-day massive clusters (1O 15 M 
by z ~ 0.5). We present an extensive spectroscopic campaign carried out on the supercluster field 
resulting, in conjunction with previously published data, in 1961 high confidence galaxy redshifts. We 
find 302 structure members spanning three distinct redshift walls separated from one another by ~ 65 
Mpc (Az = 0.03). The component clusters have spectroscopic redshifts of 0.901, 0.905 and 0.905. 
The velocity dispersions are consistent with those predicted from X-ray data, giving estimated cluster 
masses of ~ 10 14 ' 5 — 10 14 ' 9 M Q . The Dressier- Shectman test finds evidence of substructure in the 
supercluster field and a friends-of-friends analysis identified 5 groups in the supercluster, including 
a filamentary struct ure stretching between two cluster cores previously identified in the infrared by 
iCoppin et al.l (|2012l ). The galaxy colors further show this filamentary structure to be a unique region 
of activity within the supercluster, comprised mainly of blue galaxies compared to the ~43-77% red- 
sequence galaxies present in the other groups and cluster cores. Richness estimates from stacked 
luminosity function fits results in average group mass estimates consistent with ~ 10 13 M Q halos. 
Currently, 22% of our confirmed members reside in > 10 13 Mq groups/clusters destined to merge 
onto the most massive cluster, in agreement with the massive halo galaxy fractions important in 
cluster galaxy pre-processing in N-body simulation merger tree studies. 

Subject headings: Galaxies: clusters: individual (RCS 231953+0038.0, RCS 232002+0033.4, 
RCS231948+0030.1) Galaxies: high-redshift Spectroscopy: galaxies 



1. INTRODUCTION 

In a universe dominated by hierarchical structure for- 
mation, galaxies accrete along high density filaments and 
sheets to form groups and clusters. Superclusters repre- 
sent the largest density enhancements in this cosmic web 
and can reach scales of ~ 100 — 200 h' 1 Mpc, with close, 
aligned clusters ha ving the densest filamentary cluster- 
clust e r bridges (e.g. iBqnd et al.lll996t iBahcall fc Soneiral 
[T98llWrav et al.H2006lL ~ 

As superclusters represent the largest component of 
the large-scale structure in the univers e they can be 
used to test cos mological models (e.g. IBahcalll 119881 : 
iKolokotronis et al J 12003 lEinasto et al.ll201lD . They also 
offer a range of spatial and dynamical sub-environments, 
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from infalling galaxies, groups and filaments to the high- 
density cluster cores, and are therefore ideal laboratories 
in which to study the effects of the environment on galaxy 
evolution. In the nearby unive rse superclusters have 
long been accessible for study Ce.g.rGregorv fc Thompson! 
1978; IGregorv et aI1H98lt lQori l983) and catalogs of su- 
perclusters in the low-redshift universe have been com- 
piled through the large spectroscopic cam paigns done by 
the S l oan Digital Sky Survey (SDSS) (e.g. lEinasto et al.l 
I2002t iLiivamagi etaLl 120121 ) and the 2 degree Field 
Gala xy Redshift Survey (2dFGRS) (e.g. lEinasto et all 
120071) . 

The study of superclusters in the high redshift uni- 
verse, however, remains elusive. In N-body simulations 
of hierarchical structure formation, the density of super- 
clusters in the universe decreases with redshift, with few 
rich, compact su perclusters found at high redshift (e.g. 
IWrav et all 120061). Though t he number of superclusters 
found bv iWrav et al.l (|2006l ) decreased from z ~ to 
z ~ 1 by only a factor of a few, their mean cluster sep- 
aration increases, making them difficult to find observa- 
tionally since typical high redshift survey sizes are not 
large enough to recognize these immense structures and 
superclusters on smaller scales become more rare. 

To date, only a handful of spectroscopically confirmed 
superclusters have been found at high redshift. The 
CL0016 supercluster at z = 0.55 extends over 20 Mpc 
and includes three clusters along with compan ion groups 
and filamentary structure (|Tanaka et al.l l2009) . At z ~ 1 
three superc l usters of varying scales have been reported. 
iLubin et al.l (|2000l ) identified a supercluster containing 
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ten spectroscopically confirmed cluster and group candi- 
dates within a radial length of ~ 100 hjQ Mpc and trans- 
verse length of ~ 13 hm M pc at z = 0.9 (iGal fc Lubinl 
[200l IGal et"all[2008h . The ISwinbank et al.l (I200l ~su^ 
percluster at z = 0.89, found in the UK Deep In- 
frared Sky Survey (UKIDSS) Deep eXtragalactic Survey 
(DXS), contains five clusters (~ lO 13 5 ^ 14 M Q ) spread 
over 30 Mpc on the sky. The RCS 2319+ 00 supercluster 
was fir st spectroscopically confirmed by Gilbank et al. 
(|2008aT> and contains three massive clusters at z = 0.9 
with X-ray masse s of M x ,tot ~ 4.7 - 6.4 x 10 14 M Q 
(|Hicks et al.l [20 08). The highest redshift supercluster 
confirmed to date is the Lynx Supercluster of two X- 
ray emitting clusters at z = 1. 26 and z = 1.27, first 
discovered bv lRosati et al.l ([1999D . with seven additional 
candidate group s and clusters ident ified through photo- 
metric redshifts (|Nakata et al.ll2005l) . 

In this paper we present the results of a large spectro- 
scopic campaign on the RCS 2319+00 supercluster field 
at z ~ 0.9. These data will form the basis of a series 
of upcoming papers studying the effects of the high red- 
shift supercluster environment on its member galaxies. 
We aim to map out the substructure of the supercluster 
to provide a spectroscopically confirmed catalog of mem- 
bers in the regions of various densities in the cluster field 
for further studies on the galaxy populations. The paper 
is organized as follows. In fJ2] we outline the currently 
published information available on the RCS 2319+00 su- 
percluster. Section [3] details the extensive spectroscopic 
campaigns that have been completed on the superclus- 
ter while 2] summarizes the output of the redshift cata- 
log. In ^5] we attempt to quantify the cluster properties 
including the redshifts, velocity dispersions, virial radii 
and masses of the three main component clusters. We 
further trace the substructure of the supercluster in <JBJ 
The recovered structure is used to investigate the color- 
magnitude relation and galaxy color distribution in dif- 
ferent density environments in SjT] In f}8] we attempt to 
place constraints on the masses of the potential infalling 
groups identified in $6] and put them in context with cur- 
rent theories of group infall in clusters. Finally, a sum- 
mary is offered in §§\ Throughout this paper we assume 
cosmological para meters of 17a = 0.73, f2 m = 0.27, and 
H = 71 km s" 1 (jLarson et alJl2011l ). 

2. THE RCS 2319+00 SUPERCLUSTER 

RCS 231953+0038.1 (hereafter referred to as cluster 
A) at z ~ 0.9 was iden tified in the first Red- Sequence 
Cluster Survey (RCS-1; IGladders fc ^l200l a s a re- 
markable strong-lensing cluster (|Gladders etaLI I2003T) . 
Subsequent investigation of the RCS-1 catalog identified 
two potential companion clusters (RCS 232003+0033.5 
- cluster B and RCS 231946+0030.6 - cluster C) whose 
red-sequence peaks fall at a similar photometric redshift 
to cluster A. The density profiles of the red-sequence 
galaxies indicate that the three clusters are separated 
by < 3 Mpc in the plane of the sky with each cluster core 
seemingly aligned towards its closest neighbor, further 
substantiating the claim that these clust ers are bound to 
each other in a supercluster (as found by[Binggeli (jl982| ) 
in a study of Abell clusters). 

Chandra X-ray imaging confirmed the presence of 
three X-ray luminous cluster cores with Lx ~ 3.6 — 



7.6 x 10 44 erg s" 1 and M x ,tot ~ 4.7 - 6.4 x 10 14 M© 
(|Hicks et al.ll2008D . Spectroscopic confirmation that the 
three clusters all lie at z ~ 0.9 was first s hown in 
iGilbank et al.l ()2008al ). Using the meth ods of iSarazin l 
(12002ft to estimate the merger rate, Gil bank et al.1 
(|2008aj ) postulate that the three clusters will merge to 
form a Coma-like cluster (~ 10 15 Mq) by z ~ 0.5. 

As the two other confirmed z ~ 0.9 s upercl usters 
of ISwinbank etal] (|2007f ) and iLubin et al.l (j2000f ) have 
lower mass cluster components that are spread over 
larger scales, they are therefore les s likely to merge into a s 
massive a cluster as RCS 2319+00 ()Gilbank et al j2008aft . 
Thus RCS 2319+00 is a rare example of a massive cluster 
caught in the process of forming through the merging of 
multiple cluster and group components and is an ideal 
laboratory in which to study the role of the environment 
and hierarchical processes in galaxy evolution. 

Extensive follow-up is underway on the RCS 2319+00 
supercluster to study the effect of this unique environ- 
ment on its member galaxies. Cluster A has already 
been stu died in surveys of weak-lensing in massive galaxy 
clusters (pee et al.ll201lD . the Su nyaev-Zel'doyich (S Z) ef- 
fect of strong lensing clusters ([Gralla et al.l 1201 ID and 
the s ubmillimetre sour ce counts of high-redshift clus- 
ters (|Noble et al.l 120121 ) . The X-ray properties of the 
three component clusters of the RCS 2319+00 super - 
cluster have been i nvestigated by IHicks et al.l (2008). 
iCoppin et al.l (|2012l ) discovered an infrared-bright fila- 
ment of galaxies connecting component clusters A and 
B with a total SFR~ 900 M Q yr" 1 . This data paper is 
the first in a series to further investigate the hierarchical 
build-up of large scale structure at z ~ 1 through the 
RCS 2319+00 supercluster environment. 

3. THE SPECTROSCOPIC DATA 

The spectroscopic catalog emanates from six different 
observing runs with five different Multi-Object Spectro- 
graphs (MOS): Magellan-IMACS, VLT-VIMOS, Subaru- 
FMOS, Gemini-GMOS, and VLT-FORS2. TableQ]sum- 
marizes the central coordinates of each of the spectro- 
scopic masks, the dates or program ID of the observa- 
tions, the total integration time, and the number of tar- 
get slits per mask. 

The observation setup and data reduction techniques 
used for each instrument are detailed below. Figure Q] 
shows the coverage of the IMACS, VIMOS and FMOS 
masks over the supercluster field. The cluster positions 
are identified by their X-ray contours (see Figure 1 of 
Hicks et al. 2008 for a more detailed view of the X-ray 
flux distribution). The GMOS and FORS2 fields of view 
are smaller and fall within the larger instrument masks, 
centered on cluster A. 

3.1. IMACS Data 

Wide field MOS was taken with the Inamori Magel- 
lan Areal Camera and Spectrograph (IMACS) on the 
Baade 6.5 meter Magellan telescope using the Univer- 
sity of Toronto's share of Magellan time. Three pointing 
positions of the 27' diameter IMACS field of view were 
used. The data was taken as part of a large optical spec- 
troscopic survey covering 42 RCS1 and 1 RCS2 clusters 
over eight observing runs from October 2004 to October 
2007 (R. Yan et al. 2012 in preparation). Three masks 
were centered on the RCS 2319+00 supercluster field 
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Table 1 

Spectroscopic mask observations on the RCS 2319+00 supercluster field. 



Instrument 


Mask 




RA 


DEC 


Date/ID 


Total lnt time (s) 


# slits 


IMACS 


231953+A 


23 


19 


53.00 


+00 


•is 


00.00 


October 29, 2006 


7200 


722 


IMACS 


231953+B 


23 


19 


53.00 


+00 


38 


00.00 


October 9, 2007 


7200 


708 


IMACS 


231953+C 


23 


19 


53.00 


+00 


38 


00.00 


October 10 2007 


7200 


644 


IMACS 


2319NE 


23 


18 


31.00 


+00 


34 


12.00 


October 9,' 2007 


7200 


751 


IMACS 


2319W 


23 


21 


02.00 


+00 


50 


30.00 


October 10, 2007 


7200 


740 


VIMOS 


RCS2319P1 


23 


19 


32.93 


+00 


•12 


26.10 


July 9, 2005 


2500 


341 


VIMOS 


RCS2319P2 


23 


19 


32.43 


+00 


33 


55.44 


August 13, 2005 


2500 


337 


VIMOS 


RCS2319P3 


23 


20 


13.13 


+00 


32 


11.90 


August 22, 2006 


1920 


355 


VIMOS 


RCS2319P4 


23 


20 


13.70 


+00 


11 


09.60 


October 16, 2006 


2880 


339 


FMOS 


RCS2319P1 


23 


19 


49.83 


+00 


31 


41.90 


October 16, 2011 


7200 


171 


FMOS 


RCS2319P2 


23 


19 


54.30 


+00 


35 


11.60 


October 16, 2011 


7200 


171 


GMOS-N 


RCS2319+0039 


23 


19 


53.40 


+00 


38 


13.80 


September 30, 2003 


10800 


27 


FORS2 


RCS231953+0038.0 


23 


19 


53.00 


+00 


38 


00.00 


70.A-0378(B) 


10800 


29 


FORS2 


R2319M01 


23 


19 


53.00 


+00 


37 


13.00 


October 30, 2011 


1800 


21 


FORS2 


R2319M01 


23 


19 


53.00 


+00 


39 


27.00 


November 1, 2011 


2700 


20 



Note. — Units ol right ascension are hours, minutes, and seconds and units of declination are degrees, arcminutes, and arcseconds. 




23 h 21 m 00 s 20 m 30 s 00 s 19 m 30 s 00 s 



Right Ascension 

Figure 1. Spectroscopic mask coverage centered on 

23 h 19 m 52 s .8 + 00°34'12" over a 35 x 35 arcminute region 
with the three comp onent cluster posit ions identified by the 
X-ray contours from Hicks et al. (2008). Rectangles denote 
the four VIMOS pointings with each color representing the 
four quadrants of a single VIMOS observation. The orange 
circles show the three IMACS positions and green circles rep- 
resent the FMOS fov. The coverage with CMOS and FORS2 
is not shown. Both instruments were centered on cluster 
A, with CMOS targets covering ~1 square arcminute and 
FORS2 targets over a ~6 square arcminute region. 

(masks 231953+A, 231953+B, and 231953+C). Two ad- 
ditional masks covering the nearby RCS 231831+0034.2 
(2319NE) and RCS 232102+0050.5 (2319W) RCS 1 clus- 
ter candidates overlap the supercluster field region. Each 
mask was observed in 4 x 1800s exposures with ~ 650- 
750 slits per mask. The gl50 grism was used with the f/2 
camera and 1 arcsecond slit widths, yielding a spectral 
resolution of R=555. A custom band-limiting filter was 
used to restrict wavelength coverage to 6050-8580 A en- 
compassing the [OII]A3727 to G-band (~4304 A) spectral 
features and slightly beyond at the supercluster redshift. 
In total, there are 2138 unique IMACS spectroscopic tar- 
gets in a 35 square arcminute region around the super- 



cluster, with duplicates used to check the redshift ac- 
curacy of the IMACS sample (see R. Yan et al. 2012, in 
preparation) . 

The full details of the IMACS reduction and red- 
shift finding are discussed further in the R. Yan et al. (in 
preparation) paper detailing the large IMACS RCS sur- 
vey. Briefly, the IMACS spectroscopy was reduced us- 
ing a modifie d version of the D EEP2 DEIMOS reduc- 
tion pipeline (jDavis et all I2002D . The redshift finding 
was done using the IDL-based zspec program, which uses 
chi-squared fitting to template spectra. It was adapted 
from the DEEP2 software for use with the RCS IMACS 
spectroscopy sample by R. Yan. All redshift assignments 
were checked by eye and assigned a quality flag: Q > 3 
for a secure redshift, Q = 2 for a potential redshift, Q 
= 1 for bad spectra, and Q = -1 for stars. Using high 
quality spectra of repeated targets, the redshift error for 
the IMACS sample was calculated to be Az ~ 0.000388, 
or ^60 km s _1 at z ~ 0.9. 

3.2. VIMOS Data 

Optical spectroscopy was take n with the Visible Mu lti- 
Object Spectrograph ( VIMOS: Ee Fevre et al1l2003D on 
the 8.2 meter Unit Telescope 3 of the Very Large Tele- 
scope (VLT UT3). VIMOS contains four detectors, each 
with a 7x8 arcminute fov, separated by 2 arcminute 
gaps between quadrants. Four VIMOS pointings were 
designed to cover a ~25 square arcminute region centered 
on RCS 2319+00 cluster A (see Figure [J). At the time 
of the spectroscopic proposal, clusters B and C had not 
yet been identified. No color or magnitude selection was 
applied to the mask design; however, the highest prior- 
ity wa s given to sources with radio and submm counter- 
parts ([Noble et al.ll2012T ). The low resolution red grism 
(LR_red), with a wavelength rang e of 5500-9500 A, was 
used to cover the [OII]A3727 to [OIII]A5007 spectral fea- 
tures at z ~ 0.9. This set-up resulted in a spectroscopic 
resolution of R=210 with a dispersion of 7.3 A pixel -1 . 

The VIMOS spectroscopy was run through the stan- 
dard VIMOS data reduction pipeline. We decided to re- 
extract the ID spectra from the 2D sky-subtracted, flat- 
fielded, wavelength-calibrated spectra in order to (a) ver- 
ify that the correct primary object spectra was extracted 
when multiple sources appeared in a single slit, (b) im- 
prove the trace of faint continuum sources and (c) allow 
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for better removal of contamination from cosmic rays, 
th , -I s * and 2 nd order contamination and fringing in the 
spectra. The IRAF noao package was used to trace and 
extract the ID spectra, with additional background sub- 
traction performed and bad pixels replaced using model 
points. Spectra of separate observations were then coad- 
ded with imarith and checked for any residual bad pixels, 
which were removed manually. 

To prepare the spectra for redshift finding the contin- 
uum was fitted and subtracted using fitld and stacked 
with a weighted spectrum to remove the A-band absorp- 
tion line. Redshifts were found by first running the final 
ID spect ra through a cr oss-correlation program adapted 
for IDL (|Blindertl 12001 ) from the CNOC1 techniqu e us- 
ing five template galaxy spectra (|Yee et al.l I1996T) and 
then interactively checking and, if necessary, re-fitting 
the redshifts using an IDL-based redshift program writ- 
ten for use with the output of the cross-correlation pro- 
gram. The redshifts were assigned the same quality flags 
as the IMACS data. The VIMOS redshifts carry errors of 
Az ~ 0.002 (—300 km s -1 ) based on the wavelength cali- 
bration and dispersion of the data, the widths of the cor- 
relation peaks and matched targets between the IMACS 
and VIMOS samples. 

3.3. FMOS Data 

Near-Infrared spectroscopy was obtained with the 
Fiber Multi Object Sp ectrograph (FMOS) on the 8.2 me- 
ter Subaru telescope (|Kimura et al.ll20"io() . FMOS con- 
tains 400 fibres with an accuracy of 0.2 arcsec RMS in 
a 30 arcminute diameter field of view. Two pointings 
were centered on the supercluster field. We used the 
Cross-Beam-Switching (CBS) mode, which uses nodding 
for better sky subtraction by allocating two fibre spines 
per target: one that is initially on target and one that 
is initially on the nearby sky. This results in only <200 
science targets per pointing but allows each target to 
be observed for 100% of the exposure time and provides 
more accurate sky subtraction. We used the Low Res- 
olution mode (R=600, AA ~ 20 A) with a wavelength 
range of 0.9-1.8/im, which allows coverage of H Q (6563 A) 
and the SI (—6716/6730 A) doublet features and beyond 
at the supercluster redshift. Spectroscopic targets were 
chosen based on our CFHT WIRCam J and Ks catalogs, 
with priority going to mid-infrared 24/im counterparts 
from our Spitzer MIPS catalog (T.M.A. Webb et al. 2012, 
submitted). Previously identified (Q>3) and suspected 
(Q=2 redshifts) luminous infrared galaxy (LIRG) super- 
cluster members from the IMACS and VIMOS redshift 
lists were targeted for further confirmation of their red- 
shifts and for use in further studies of their spectral fea- 
tures. Overall 345 targets, including 121 targets not 
previously observed by other instruments, were observed 
with FMOS. 

The FMOS spectroscopy was run through the FMOS 
preliminary data reduction and calibration package, 
FIBRE-pac: FMOS Image-Based REduction Package 
(jlwamuro et al.l [2012f ). which uses IRAF scripts and C 
programs based on CFITSIO to produce ID spectra. The 
reduction includes flat-fielding, bad pixel removal, spa- 
tial distortion correction, spectral distortion correction, 
background subtraction, bias correction, wavelength cal- 
ibration and flux calibration. The resulting spectra have 



a dispersion of 5 A pixel -1 . Final ID spectra were pre- 
pared in the same way as the VIMOS spectra for use with 
the cross-correlation and redshift confirmation programs, 
with redshifts having the same quality flags and errors as 
the previous data. The redshift errors were again taken 
to be Az ~ 0.002 to account for the wavelength disper- 
sion and correlation peak widths. For any target with 
a previous, low confidence spectra from IMACS or VI- 
MOS, the estimated FMOS redshift was checked with 
the lower wavelength range spectra in an effort to con- 
firm the FMOS redshift and/or increase the confidence 
of the previous redshift. As FMOS is still in shared risk 
mode and the FIBRE-pac reduction package is still in 
development, the data were noisier than anticipated, de- 
spite photometric observing conditions, with only —30% 
redshift success. 

3.4. GMOS North Data 

Spectra from the 5.5 arcminute field of view Gem- 
ini Multi-Object Spectrograph (GMOS) on the 8.1 me- 
ter Gemini North Telescope was obtained in an effort 
to identify the redshift for the strong lensing arc in 
cluster A (jGilbank et al.l I2008al) . The data was ob- 
served in the band-shuffle nod-and-shuffle mode, result- 
ing in 27 targets in a — 1 square arcminute region. The 
R150_G5306 red grism was used with a band-limiting 
filter (GG455 G0305) and 2x2 detector binning for a 
wavelength range of -4500-10000 A, allowing coverage 
of [OII]A3727 to [OIII]A5007 at z - 0.9 and giving a 
resolution of 11.48 A and a dispersion of 3.564 A pixel -1 . 

The GMOS spectroscopy was reduced in IRAF with 
the help of the gemini package tasks to bias-subtract, flat- 
field, and wavelength-calibrate the data. The positive 
and negative spectra for each object from the nod-and- 
shuffle observations were coadded and ID spectra were 
extracted using the IRAF a pa 1 1 task with residual bad 
pixels removed manually. The spectra were prepared and 
redshifts found as in the VIMOS and FMOS data, with 
the same errors and quality flags assigned. 

3.5. FORS2 Data 

Spectroscopy taken with the FOcal Reducer and 
low dispersion Spectrograph (FORS2) of the Very 
Large Telescope was obtained from two separate ob- 
servation runs. In the first observation run (mask 
RCS231953+0038.0) spectra of 29 target galaxies was 
taken as part of the larger spectroscopic follow-up of 
RCS-1 clusters. FORS2 was used in MXU multi-object 
spectroscopy mode with the 3001 grism in conjunction 
with the OG590 filter to give a dispersion of 108 A mm -1 
and wavelength coverage from 6000-11000 A. The target 
galaxies were selected primarily using a fixed magnitude 
limit. The observations totalled three hours per mask, 
including overheads, split into individual observations to 
improve cosmic ray rejection. Redshifts were found us- 
ing cross-correlation with the IRAF rvsao package and 
interactively checked and, where necessary, refitted us- 
ing the emsao task. Quality flags have been adjusted to 
be consistent with our larger redshift catalog. 

The second FORS2 data set consisted of a total of 43 
slits over two masks (R2319M01 and R2319M02) taken 
as a follow up of Hers chel selected sourc es from the Her- 
schel Lensing Survey (|Egami et al.ll2010[ ). The data were 
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taken using the 300V grism combined with the GG435 
filter to give a wavelength coverage of 445-865 A with 
a 112A mm dispersion. Rcdshifts were found manu- 
ally using the peak of the [OII]A3727 emission line. The 
targets in this run comprised mainly background galax- 
ies lensed by cluster A of the RCS2319+00 supercluster, 
however six fell within the supercluster redshift range, 
with three having no previous spectroscopy. 

4. THE RCS 2319+00 REDSHIFT CATALOG 
4.1. The combined spectroscopic data set 

In total, our spectroscopic campaign in the 35 x 35 ar- 
cminute field around the RCS 2319+00 supercluster con- 
tains 3305 unique spectroscopic targets. Where there are 
duplicate targets present between masks or different in- 
struments, the redshift with the higher quality factor is 
used. If the quality factor is the same between spectro- 
graphs, then priority is given to IMACS redshifts since 
it is the instrument with the smallest error. In total, 
there are 66 duplicate targets between the IMACS and 
VIMOS samples that have confident redshifts with both 
instruments. Of these, 9% had redshifts that did not 
fall within 3(7 of each other. This implies that, in the 
worst case scenario, 8 of the VIMOS structure members 
listed in Table [2] have been mis-identified. The average 
redshift error on the IMACS/VIMOS matched galaxies 
is Az ~ 0.002. The redshifts of overlapping targets be- 
tween the combined IMACS/VIMOS spectroscopic cat- 
alog and the FMOS data were found with the aid of 
the previous spectra and consequently they all match 
within errors. The redshifts of the duplicate sources in 
the FORS2 spectroscopy matched to the combined red- 
shift list to within Az ~ 0.001 and the two duplicates 
from the GMOS sample matched to within Az ~ 0.003. 

The full redshift distribution of the RCS 2319+00 su- 
percluster region is shown in the top panel of Figure [21 
Our catalog contains 2000 good confidence redshifts, 
with 1961 galaxy spectra and 39 stars. There is a clear 
peak at the supercluster redshift of z ~ 0.9 along with 
potential foreground peaks at z ~ 0.28, 0.45, 0.53 and 
0.76. The photometric redshifts of the RCS 2319+00 
clusters, as identified by overdensities of red galaxies in 
this field (L.F. Barrientos et al. 2012, in preparation), are 
shown as dashed lines in the top panel of Figure[2](A & B: 
z p hot,RS = 0.984, C: z phot , RS = 0.895). The red-sequence 
redshifts were found to carry errors of Az < 0.05 over 
the majority of the RCS-1 survey range, increasing to 
Az ~ 0.1 at the highest redshifts due to larger uncer- 
tainties in the photometric data and poorer sampling of 
the 4000 A break in the R c and z' bands. 

The bottom panel of Figure [2] shows the zoomed 
in redshift histogram covering a redshift range of 
0.85 < z < 0.95. Three distinct redshift peaks are 
apparent, separated by ~ 4700 kms -1 (Az ~ 0.03) from 
one another, ~ 65 Mpc assuming Hubble flow. To deter- 
mine membership in the RCS 2319+00 supercluster field 
we use a broad membership range of 0.858 < z < 0.946 
based on the bottom panel of Figure [21 This is done to 
encompass any infalling groups, filaments or walls in red- 
shift space. This is the range over which we search for 
substructure in the supercluster and may contain isolated 
galaxies not connected to the supercluster structures (see 
Ell- 



Table 2 

Spectroscopic target numbers in the RCS 2319+00 supercluster 

field 



Spectroscopic 


Targets 


Secure 


Stars 


Galaxies 


Structure 


Source 




Rcdshifts 






Members 


IMACS 


2133 


1325 


33 


1292 


185 


VIMOS 


962 


562 


6 


556 


85 


FMOS 


181 


87 





87 


25 


GMOSN 


14 


11 





11 


4 


FORS2 


15 


15 





15 


4 


TOTAL 


3305 


2000 


39 


1961 


302 



Note. — All numbers listed are for unique spectroscopic targets 
only. Any duplicate targets have been attributed to the instrument 
with the highest redshift quality and/or the lowest redshift error. 



A different view of this redshift range is shown in Fig- 
ure[3l where each of the three panels represents the veloc- 
ity versus radius of all spectroscopically confirmed galax- 
ies from the center of each cluster relative to z — 0.9. The 
dotted lines indicate the positions of the ± 1 Mpc radii 
from the neighboring cluster centers relative to their ra- 
dial position from the center of the specified cluster. The 
radial proximity of the clusters highlights the difficulty 
in assigning cluster membership over the field to specific 
cluster cores. The three separate velocity walls can be 
seen in all three panels at z ~ 0.87, 0.90, 0.93. These 
redshift sheets may or may not belong to a single, grav- 
itationally bound structure. For the purpose of investi- 
gating substructure around the supercluster region, we 
retain all spectroscopically confirmed galaxies within the 
broader redshift range. 

Table [2] gives a summary of the number of unique spec- 
troscopic targets and the good confidence redshifts for 
each spectrograph along with the number of good confi- 
dence stars, galaxies and overall structure members. For 
the rest of this paper we use only secure spectroscopic 
rcdshifts. We find 302 members in the broad superclus- 
ter range of 0.858 < z < 0.946. 

4.2. Spectroscopic completeness and success rates 

Due to the large survey area of the supercluster field, 
targeting all sources with spectroscopy is unrealistic. 
The 5(7 magnitude limits of th e RCS survey are z' = 
24.05 mag and R c = 24.7mag (jGladders fc Yed 120051) . 
resulting in costly time overheads to obtain high S/N 
spectroscopy for the faintest sources in our catalog. From 
Figure [T] it can be seen that the spectroscopic coverage 
of the three cluster cores varies, with cluster A having 
the most mask coverage from all instruments including 
the additional GMOS and FORS2 data. Cluster B has 
the least amount of spectroscopic coverage, with its core 
falling outside or along the edges of the VIMOS mask 
coverage. 

The completeness of the spectroscopic target coverage 
over the entire optical catalog for magnitude bins of 0.5 
in z' magnitude over the 35 x 35 arcminute area is shown 
in the left panel of Figure [H The solid line shows the 
target completeness for all spectroscopic targets, regard- 
less of redshift success, while the dashed line shows only 
those targets that yielded a secure redshift. Despite our 
3305 independent spectroscopic targets, our highest tar- 
get completeness, for 20-21 magnitude sources, is only 
^37%, decreasing to ~30% coverage for secure redshifts. 

Spectroscopic selection functions and weights were 
calculated for each source with a spectroscopic red- 
shift following the method used in the two Canadian 
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Figure 2. Top: Spectroscopic redshift distribution in the 35x35 arcminute region around the RCS 2319+00 supercluster field 
with a bin size of z = 0.01 (~ 1500 km s _1 at z ~ 0.9). The filled histo gram represents good confidence redshifts while the 
white histogram shows all redshift results. Potential foreground structures are visible in the histogram, as well as a clear peak at 
z ~ 0.9. The dashed lines indicate the photometric redshifts for the three component clusters of RCS 2319+00 (A & B: z p hot,RS 
= 0.984, C: z p hot,RS = 0.895; L.F. Barrientos et al. 2012, in preparation). The la error on the RCS-1 photometric redshifts at 
z ~ 0.9 is Az ~ 0.1 and is shown by the horizontal error bars. Bottom: Zoomed in redshift histogram around the supercluster 
redshift of z ~ 0.9 with a bin size of z = 0.002 (~ 300 km s _1 at z ~ 0.9), consistent with the largest redshift error in our 
catalog. 



Network for Obser v ational Cos mology (CNOC) sur- 
veys (see lYee et all (|1996l [2000) for a detailed expla- 
nation). Briefly, the magnitude selection function (S m ) 
of each spectroscopically confirmed supercluster member 
is found by taking all galaxies within the photometric 
catalog in a magnitude bin centered around the desired 
galaxy and calculating the ratio of galaxies in that bin 
having a measured redshift to those without a redshift. 
The geometric distribution selection function {S xy ), used 
as a correction to the primary magnitude selection func- 
tion, was also computed. The geometric selection func- 
tion corrects for different target densities over the larger 
field and uses only sources within 2 arcmins of the de- 
sired galaxy, within ±0.5mag bins to compute a local 
magnitude selection function. This is then divided by 
the previously found S m . The spectroscopic weight of 
each galaxy is the inverse of its total selection function 

("5* = S m S X y). 



The spectroscopic redshift success rates for the tar- 
geted sources is shown in the middle panel of Figure [4] 
and represents all instruments combined, with ^80% of 
the brighter targets yielding confident redshifts and suc- 
cess rates dropping to only ~30% at our 5<r depth of z' 
= 24.05mag. The S/N, wavelength dispersion and re- 
duction methods all effect the quality of the final spectra 
and therefore the redshift success rate. Most spectra 
were not taken under photometric conditions and the 
S/N varies between different spectra from the same in- 
strument based predominantly on the weather conditions 
throughout the observations. The success rate of bright 
(z' < 22mag) and faint (z' > 22 mag) sources, based 
on R c -z' color is also shown (Figure [4j right). For bright 
sources, the redshift success is fairly consistent at 70- 
80% for all colors, with the lowest and highest color bins 
affected by small numbers. 
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Radius from Cluster Center (Mpc) 
Figure 3. Velocity versus radius of all spectroscopically con- 
firmed galaxies in the RCS 2319+00 supercluster field within 
the redshift range 0.848 < z < 0.946. In each panel the 
galaxies' radial positions are plotted relative to the specified 
cluster center. All velocities are plotted relative to an arbi- 
trary z = 0.9. The dotted lines denote the ± 1 Mpc radii for 
the neighboring cluster centers relative to their radial posi- 
tion from the center of the specified cluster to demonstrate 
the difficulty in assigning cluster membership to specific clus- 
ter cores. 

5. CLUSTER PROPERTIES 
5.1. Determining the proper cluster centers 

Figure [5] shows the distribution of galaxies within the 
RCS 2319+00 supercluster field over the 0.858 < z < 
0.946 range discussed above, plotted over the density 
profiles of the red-sequence galaxies map at a photo- 
metric redshift of z ~ 0.9 (L.F. Barrientos et al. 2012, 
in preparation). The symbol size and color varies with 
the redshift of the galaxies: larger, darker symbols are 
blue shifted relative to the mean redshift; smaller, lighter 
symbols are red shifted. The X -ray contours seen in Fig- 
ured] ffrom lHicks et all 120081 ) are not plotted but their 
peaks are shown with large green X marks and show an 
offset to the red-sequence peaks (see Table [3] for coordi- 
nates and offset distances) . These offsets are most likely 
due to the contamination from the lower redshift wall 
seen in Figures [2] and [3] causing the cluster peak posi- 
tions in the red-sequence density profile, which spans a 
broad redshift bin, to shift. The RCS 2319+00 superclus- 
ter is already visible in the red-sequence density profile 
maps with central redshifts as low as z ~ 0.758. As our 
spectroscopic coverage is not complete, the galaxies in 
this lower redshift peak, whose spatial distribution can 
be seen in Figure [5] represent a lower limit on the red- 
sequence density profile contamination. 

We look to the brightest cluster galaxy (BCG) in each 
cluster core to aid in isolating the correct cluster cen- 
ter. The BCGs were found by first isolating all galaxies 
less than 0.5 Mpc from the red-sequence cluster centers 
in the photometric catalogue. From the reduced source 



lists, the brightest z'-band galaxy within ±0.2 of the 
R c — z' red-sequence line for each cluster was selected as 
the BCG. The identified BCGs for the three clusters are 
all found <6.5 arcseconds from the cluster X-ray peaks, 
closer than the 13.5-32.5 arcseconds between the X-ray 
and red-sequence cluster center positions. We therefore 
use the X-ray centers in calculating the dynamical prop- 
erties of the clusters. 

5.2. Determining the cluster core members 

To calculate the velocity dispersions and spectroscopic 
redshifts of the component clusters f §5.3|) we use all con- 
fident redshifts within radii of 0.5 and 1.0 Mpc of the 
X-ray peak positions in each cluster. Due to the close 
separation distance between the three component clus- 
ter centers, a 1.0 Mpc radius is the largest used to avoid 
double counting of galaxies between component clus- 
ters. The green dotted circles in Figure [5] shows the 
1.0 Mpc radii from the X-ray centers. An initial cut of 
0.85 < z < 0.95 around the supercluster redshift peak 
was used and the redshift histograms for the three cluster 
cores are shown in Figure [6] All secure redshifts within 
1.0 Mpc of the center positions are shown as hatched his- 
tograms, with the solid histograms including only those 
galaxies within 0.5 Mpc of the centers. The histograms 
were examined for obvious outliers not consistent with 
a Gaussian curve that could falsely broaden the velocity 
distribution. 

For cluster A, no sources were removed as the redshift 
distribution, while broad, does not show any obvious ir- 
regularity and an iterative interloper clipping method is 
used in H5. 31 when determining the cluster redshifts and 
velocity dispersions. We again look to the BCG of clus- 
ter A for comparison to ensure that the redshift distri- 
bution in the region around the cluster core is tracing 
the proper cluster redshift. The BCG is found to have a 
confident spectroscopic redshift of 0.901, consistent with 
the approximate peak of the redshift distribution in the 
left-hand panel of Figure [6] 

Clusters B and C seem to show double peaks in their 
redshift distributions, with lower redshift peaks at z ~ 
0.868 and z ~ 0.877 respectively, and higher redshift 
peaks consistent with the overall supercluster redshift 
peak at z ~ 0.9 visible in all three cluster cores. This 
could possibly be due to under sampling of the cluster 
cores, especially for B, which falls at the edges or in the 
gaps of the VIMOS spectroscopy masks (FigureQ]). How- 
ever, if this were correct and the true redshift of the 
cores falls between the two peaks, the resulting velocity 
dispersions would lead to cluster virial masses much too 
large to be consis tent with the measured X-ray masses 
(jHicks et al.l 120081 ). Another, more probable possibility 
is that the lower redshift peaks in clusters B and C be- 
long to filaments and groups of the potential foreground 
structure seen in Figure [3] We investigate this further by 
looking at the spatial distribution of all members having 
redshifts consistent with the separate peaks in the 1.0 
Mpc regions around clusters B and C. As well, we use the 
spectroscopic redshifts of the BCG in each cluster core to 
corroborate our choice of approximate cluster redshift. 

For cluster C, the redshift peak of 0.874 < z < 0.882 
(Figure [JJ center) shows that the members of this peak 
are inconsistent with the red-sequence significance con- 
tours of the cluster core. However, the members of the 
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Figure 4. Left: Spectroscopic target coverage completeness for the 35x35 arcminute catalog region. The solid line shows 
spectroscopic targets representing all slits, including those that did not yield reliable redshifts, compared to the full optical 
catalog with 5a depth of z' = 24.05 mag. The dotted line shows the coverage completeness when only the secure redshifts are 
taken into account. Middle: Spectroscopic redshift success for all targets based on z' magnitude. Right: Redshift success based 
on R c -z' color for bright (z' < 22 mag) and faint (z' > 22 mag) targets. 



Table 3 

Cluster center positions 

Rcd-scqucncc (RS) center X-R.ay center BCG position R.S/X-ray Offset X-ray/BCG Offset 

(arcsec) (arcscc) 

A 23 19 52.8 +00 38 03.5 23 19 53.2 +00 38 12.5 23 19 53.42 +00 38 13.65 TO 3~5 

B 23 20 02.5 +00 33 29.6 23 20 02.1 +00 32 57.6 23 20 02.33 +00 33 03.10 32.5 6.5 

C 23 19 46.1 +00 30 37.5 23 19 46.8 +00 30 14.3 23 19 46.60 +00 30 09.54 25.5 5.6 



higher redshift peak of 0.894 < z < 0.914 (Figure 
right), present in all three cluster cores, are consistent 
with the red-sequence density contours. The BCG in 
cluster C also validates the choice of the higher redshift 
peak as the correct one, with a spectroscopic redshift of 
0.899. 

The case of cluster B is more difficult to disentangle. 
While there are similar numbers of galaxies in the two 
redshift peaks, neither the higher, broader peak nor the 
lower, sharper peak at 0.866 < z < 0.876 (Figure [71 left) 
seem to have members significantly more clustered within 
the red-sequence over-density contours. Again, this is 
most likely due to the spectroscopic undersampling of 
cluster B compared to the other cluster cores. We rely on 
the cluster BCG, with a spectroscopic redshift of 0.903, 
to isolate the correct redshift peak for this cluster as the 
higher of the two at z ~ 0.9, consistent with both clusters 
A and B. 

We manually remove the lower redshift peaks from the 
histograms of clusters B and C so as not bias the cluster 
dynamics determinations. 

5.3. Cluster Redshifts and Velocity Dispersions 

The reduced galaxy redshift lists from the cluster core 
histograms, with the lower redshift peaks in clusters B 
and C removed as described, were used to determine the 
cluster redshifts and velocity dispersions from within the 
0.5 and 1.0 Mpc radii. W e follow the method described 
bv lPostman et ail (11998D . using the biweight mean esti- 
mate ([Beers et alj|1990l) to iteratively calculate the mean 
spectroscopic redshift and velocity dispersion of the clus- 
ter cores. At each iteration, all galaxy rest frame rela- 
tive velocities to the biweight mean cluster redshift are 



calculated as /S.v = c(z — z)/(l + z). Any galaxy 
falling outside of the mean dispersion by more than 3cr 
or 3500 km s _1 is removed from the list and the biweight 
mean is recalculated. The biweight mean is used as it is 
less susceptible to outlying interlopers that may other- 
wise disproportionally shift the value of the mean. 

The results of the mean redshift and velocity disper- 
sion measurements within the two different radii from 
the cluster centers are listed in Table 01 The larger num- 
ber of cluster galaxies used for the dynamical property 
measurements within 1 Mpc of the cluster centers result 
in smaller errors and more reliable values. For cluster B, 
there are not enough galaxies within a 0.5 Mpc radius of 
the center coordinate to calculate the dynamical proper- 
ties. We therefore use the values derived from the 1 Mpc 
radii. Cluster B again has the fewest galaxies within 
this radius with only 9 galaxies remaining after the iter- 
ative interloper clipping, making the velocity dispersion 
unreliable, as seen by the large error. The 15 galaxies re- 
maining within IMpc of cluster C offer a better estimate 
of the properties but would still benefit from additional 
spectrocopic members for a more dependable measure- 
ment. Cluster A has 23 galaxies remaining after inter- 
loper clipping, though there remains a significant error of 
~20% due to the broad distribution. This may indicate 
that we are not sampling a virialized population in this 
cluster and are instead witnessing a cluster merger. We 
discuss this possibility below when looking for substruc- 
ture both in the multiwavelength and spectroscopic data 
(see sections 15.51 and [6]) . 

5.4. Virial Radius and Mass 



The RCS231953 Supexcluster 



9 




Arc Minutes 



Figure 5. Redshift distribution of RCS 2319+00 supercluster members (0.858 < z < 0.946) centered on 23 h 19 m 52 s .8 + 
00° 34' 12" and plotted over the significance contours of the red-sequence galaxy density profile at a photometric redshift of 
z ~ 0.9. Spectroscopically confirmed members are plotted as circles, the size and color varying with redshift, with circle size 
decreasing with redshift and the colors varying from black through to yellow with increasing redshift. The dashed green circles 
denote the 1.0 Mpc radii around each clusters X-ray center, indicated with large green X marks and seen to be offset from the 
red-sequence significance peaks. 



Table 4 

Spectroscopic redshifts and velocity dispersions within two different radii of the cluster centers. 







0.5 Mpc 






1.0 Mpc 






N 


Zspec 


a 

(km s" 1 ) 


N 


Zspec 


(km s _1 ) 


A 


17 


0.900 ± 0.008 


1253 ± 306 


23 


0.901 ± 0.008 


1202 ± 233 


B 








9 


0.905 ± 0.007 


1154 ± 367 


C 


11 


0.906 ± 0.005 


759 ± 203 


15 


0.905 ± 0.005 


714 ± 180 



Using the measured velocity dispersions we can esti- 
mate the dynamical mass and radius of the component 
cluster cores. The cluster virial mass can be approx- 
imated by calculating the mass at which the overden- 
sity of the cluster is 200 times the critical density (p c ). 
We use the method of iCarlberg et al.l (|1997t ) to calculate 



M200 from the velocity dispersion using, 

3^3 cr 3 



Mo. 



H{z)G' 



(1) 



where H(z) = H ^/D,^ + Qm (1 + z) 3 and a is the one 
dimensional velocity dispersion. The virial radius is then 
approximated from M200 = (4/3)7rr| 00 x 200p c , which 
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Figure 6. Histograms of the spectroscopic redshifts within a 1.0 Mpc radius of each cluster's X-ray peak. Solid histograms 
show confident redshift within 0.5 Mpc of the cluster centers and hatched histograms extend the radius to 1.0 Mpc. 
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Figure 7. Distribution of galaxies in the three component cluster cores within three redshift slices corresponding to the lower 
redshift peaks in the histograms of Cluster B (left) and C (middle) and the higher redshift peak present in all three clusters 
(right). The dashed circles denote the 1 Mpc radii from the X-ray peaks. 



can be simplified as, 



^200 



73 



10 H(z)' 



(2) 



The calculated r2oo and M200 values are found in Ta- 
bleU 

Another method of determining the mass of clusters 
within r 2 oo fr om their veloc i ty dis persions is to employ 
the method of lEvrard et al.l (|2008f) . who used hierarchi- 
cal clustering simulations to determine a virial scaling 
relation for massive dark matter halos of 



Mo 



10 15 M P) 



h(z) 



/ G_DM_ 

Cl5 



1/a 



(3) 



where the logarithmic slope is a — 0.3361 ± 0.0026 and 
the normalization at a cluster mass of 10 15 Mq is 015 = 
1082.9 ± 4.0 km s _1 . The dark matter halo velocity 
dispersion is related to the velocity dispersion of galaxies 
by the velocity bias, b v = <J ga il odm- The current best 
estimate is (b v ) — 1.00 ± 0.05. We assume the non-bias 
case of b v = 1 , as done in previ ous cluster studies (e.g. 
iBrodwin etaUboiOtlSlfon et alj|2012f ). From the derived 
mass, we can compute the related r 2 oo value for each 



cluster. The results of the dynamical mass and radius 
using this method are listed in Table [5] 

The results for the cluster masses using the two meth- 
ods agree within errors; however, the masses derived us- 
ing the virial theorem are ~ 50% larger t han those found 
using the dark matter scaling relation. ICarlberg et al.l 
( 1996) found that masses measured using the virial the- 
orem are dependent on the boundary within which the 
cluster properties are determined and should be cor- 
rected with a surface pressure term, which requires a 
knowledge of the mass profile of the cluster. For cases 
where only the interior of the cluster's profile is being 
measured with no surface pressure correction, as is the 
case here, they found that the overestimation of the clus- 
ter mass should be less than 50%. 

With this in mind, and the results for the mass and 
virial radius found using the method of lEvrard et al.1 
(2008) being ~ 50% lower than those found with the 
virial theorem, we estimate that the masses of the clus- 
ters are closer to those calculated using the dark matter 
scaling relation. However, as the velocity dispersion mea- 
sured for cluster B is based on only 9 spectroscopic galaxy 
redshifts it, along with the resulting cluster properties, 
must be corroborated with the values determined for this 
cluster from additional supporting data before it can be 
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deemed as a reliable estimate. In the following section 
we compare the cluster parameters found using the spec- 
troscopic data to estimates from multiwavelength data. 

5.5. Comparison With Other Mass Estimators 

We can compare our results for the velocity disper- 
sions and mass estimates of the three clusters to those 
found by other methods where available. We us e the 
X-ray temperatures reporte d in Kicks et al.l (12008ft and 
the a - T x relation from Kue fc Wul (|2000D of a = 
10 2.49±o.O2 T £65±o.o3 tQ predict the ve l oc ity dispersions 
of the three clusters. The cluster masses ext rapolated 
direct ly from the X-ray data are also available. Uee et all 
(|2011[) calculated weak lensing values for cluster A from 
HST ACS/WFC d ata. Using a SZ scaling relation, 
iGralla et al.l (|2011l ) found an M 50 o value for cluster A 
of ~ 0.5 x the X-ray derived M 500y x value. They noted 
that the SZ derived masses are lower than those previ- 
ously published for a number of the clusters in their sam- 
ple from various multiwavelength methods and therefore 
we do not include the SZ mass estimate for cluster A in 
our discussion. Table [6] lists the known multiwavelength 
measurements of the masses and velocity dispersions of 
the three cluster cores. 

The velocity dispersions measured from the spectro- 
scopic data and those predicted from the X-ray temper- 
atures agree well within errors for clusters A and B. The 
values for cluster C are less consistent, varying by about 
1.5cr between the two values. This would suggest that 
the dynamical mass estimates based on the measured 
velocity dispersions should be reasonable, even for clus- 
ter B with it's low spectroscopic galaxy n umbers. In fact, 
the dy namical masses measured using the lCarlberg et all 
(|1997|) virial mass method, which is kn own to overes- 
timat e the cluster mass by up to 50% (jCarlberg et al.l 
1996), overestimates the X-ray extrapolated mass by only 
1.12cr and 0.74(7 for A and B and underestimates the 
mass by 0.45cr for cluster C from the dynamical mass er- 
rors only. The m asses estimated usin g the dark matter 
scaling relation of lEvrard et al.l (|2008l ) are all consistent 
with the X-ray masses, which are also based on a scaling 
relation to extrapolate t he mass based on the tempera- 
ture of the cluster cores (jHicks e t al. 2008]). The X-ray 
data therefore confirms that the dynami cal properties 
found from the spectroscopic data using the lEvrard et al.l 
(2008) dark matter scaling relation, with their large er- 
rors, are reasonable estimates for the cluster cores' mass 
and radius. For cluster A, the velocity dispersion found 
from the weak lensing analysis is lower than that found 
using the spectroscopic and X-ray data, though just con- 
sistent within errors of the two measurements. The weak 
lensing mass estimate is consistent with both dynamical 
mass estimates as well as with the X-ray mass. From the 
multiwavelength data currently available we estimate the 

Table 5 

Dynamical properties of the component clusters found with two 
different methods from spectroscopic velocity dispersions 



CarlberK et al. ( 1997) method Evrard et al. (2008) method 
»~200 M-200 ''200 A'hoo 

(Mpc) (10 14 A/ Q ) (Mpc) (10 14 Mq) 

~7L 1.82 ± 0.35 18.4 EE 10.7 1.58 EE 0.30 11.9 EE 6.9 

B 1.75 ± 0.55 16.2 ± 15.5 1.51 ± 0.48 10.5 ± 10.0 

C 1.08 ± 0.27 3.8 ± 2.9 0.94 ± 0.24 2.5 ± 1.9 



masses of the three clusters to be ~ 10 14 5 — 10 149 . 

The X-ray data also provide insight into the possibil- 
ity of merging activity in the cluster cores. There is evi- 
dence in the X-ray data that cluster B may contain sig- 
nificant substructure or have undergone a recent merger, 
as the peak of the X-ray emission is not consistent with 
the center of the exte nded X-ray emission (see Figure 
1 in Hicks et ail 120081) . The X-ray data for cluster A 
also seems to indicate the presence of some possible sub- 
structure in the radial surface brightness profile, but not 
enough to show evidence of a significant merger or verify 
if the popul ation we are samp ling is virialized or not (see 
Figure 2 in lHicks et aT1[2008l) . 

6. THE SUBSTRUCTURE OF THE RCS 2319+00 
SUPERCLUSTER 

Due to the close proximity of the three cluster cores, 
both in angular separation and redshift space, designat- 
ing galaxies as unique members of one of the particular 
cluster cores is difficult, as seen in Figure [3] We there- 
fore chose to study the cluster member properties as a 
function of their distribution within the overall structure 
and their apparent local density environment within the 
supercluster. The spectroscopically confirmed superclus- 
ter members have allowed us to begin tracing the super- 
cluster structure and identifying spatial and dynamical 
sub-environments within it. 

6.1. The Dressler-Shectman Test for Substructure 

As an initial test to identify the distribution of the 
underlying structure in the RCS 2319+00 s uperc luster 
field we performed the lDressler fc Shectmanl (|1988D (DS) 
test. The DS test looks for inconsistency between the 
locally measured velocity dispersion and velocity mean 
with the parameters measured over the whole structure 
field. Deviations from the global velocity parameters of 
the cluster indicate the likely presence of substructure. 
The DS test was used with the RCS 2319+00 data to 
begin tracing the locations of possible groups and fila- 
ments within the supercluster before attempting to as- 
sign specific galaxies to these potential groups. We follow 
iDressler fe Shectmanl ([1988) to compute, for each galaxy 
in our structure, the deviation <5 from the global mean as 

5 2 = {^^) [(viocai - vf + (a local a) 2 ], (4) 

where N nn is the number of nearest neighbors used to 
compute the local variables. We use 10 neare st neigh- 
bors, as done by IDressler fc Shectmanl ()1988f) in their 
cluster study. The A statistic is defined as the sum of all 
deviations Si and a cluster is said to have substructure if 
A./N mem i, ers > 1. To test for false positives, we shuffle 
the redshifts between galaxy positions and perform 25000 
Monte Carlo (MC) tests. If the substructure detected in 
the real data is indicative of true substructure, then the 
probability that the shuffled data has a higher A statistic 
than the real data will be low. The probability of false 
structure detection is 

N(A s huffle > A-real) , c n 

p = — — »r ' ( 5 ) 

^shuffle 

where N s h u f fie is the number of MC shuffling tests per- 
formed. 
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Dynamical properties from multiwavelength data 
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a X-ray mass estimates (M2oo,x) were extrapolated to T200 by Gilbank et al. (2008a). We use the X-ray temperatures reported in 
IHicks et al.1 12009 1 and the a-T x relation from lXue feWul <200Cft of a = lo 2 - 49±0 02 j*-68±0.03 to predict the ve i oc i ty dispersions of the 
three clusters. 

b Parameters are derived through weak lensing analysis of HST ACS/WFC data by Jee et al. (2011) and not through dynamical methods. 

The DS test was done both on the full redshift range of 
our larger structures and on the three separate redshift 
walls with ranges of ±2000 km s _1 for each to avoid 
significant overlap. For the full redshift range, we find 
A/A^ mern f, ers = 2.02, while the three redshift peaks yield 
values of A/N mem bers = 2.00,2.76 and 2.64 for the red- 
shift walls centered on z ~ 0.87, 0.90, and 0.93, indicat- 
ing the presence of substructure in each case. The Monte 
Carlo shuffled tests find P-values of < 0.006 for the full 
redshift range, and <0.05, 0.0002 and 0.003 for the three 
separate redshift walls, leading us to believe there is real 
substructure in each of the redshift peaks in our super- 
cluster field and prompting a search for the specific sub- 
environments within our spectroscopic galaxy catalogue 
in the following sections. 

The results of the test for the lower and central-rcdshift 
walls are shown in Figure [8] in DS test bubble diagrams, 
with each galaxy represented by a circle whose size scales 
with e s for that galaxy. The three cluster cores are indi- 
cated by the IMpc radius dashed circles on the each of 
the panels. Though the highest redshift slice has a high 
A/ ^members of 2.64, it has the lowest number of galax- 
ies, causing the 10 nearest neighbours to be spread over 
a large region in most cases and so we do not find any 
evidence of substructure in concentrated regions. 

In the z ~ 0.87 redshift slice, we find that both clus- 
ters B and C show deviations from the overall mean pa- 
rameters, with substructure extending from the centers 
of both cluster core positions. Deviations are found to 
extend from the center of cluster B slightly to the north- 
east as well as a more elongated deviation north-west 
towards the middle of the three clusters. These areas 
of broad redshift distribution can be seen in the top 
panel of the velocity versus radius plot in Figure [3] as 
a wide grouping of galaxies, both in radius and velocity, 
~ 1.5 — 2.5 Mpc from the center of cluster A. Whether 
the deviations that extend in the two directions outward 
from cluster B belong to the same substructure or to sep- 
arate overdensities in the supercluster will be explored in 
the following section. In cluster C, evidence of substruc- 
ture is shown to extend south west of the cluster core, 
with the largest deviations from the mean present in this 
region. 

The redshift slice around z ~ 0.9 also shows indications 
of substructure in the area of cluster C, though in a more 
concentrated region and with lower deviation than seen 
in this cluster in the lower redshift slice. There are also 
regions of potential substructure east of cluster B and 
west of clusters A and C. An interesting area to note is 
the region between clusters A and B, where there seems 
to be some evidence of substructure consistent with a 



Figure 8. DS test bubble diagrams with each galaxy in our 
structure catalogue denoted by a circle whose radius scales 
with e s . Dashed circles show the IMpc radius around each 
cluster center. Results from the DS test on the ±2000 km 
s _1 region around the lower z ~ 0.87 redshift wall (left) and 
the z ~ 0.9 redshift structure (right). 

previously confirmed infrared bright filament found in 
the Herschel 250-500/im Spect ral and Photometric Imag- 
ing REceiver (SPIRE) map bv lCoppin etaU (|20l2T) . The 
complementary detection of this structure in the IR pro- 
vides support for our claim of a substructure detection 
in this area. 

Both clusters B and C were seen to show two redshift 
peaks in Figure [5] within the central 1 Mpc around the 
clusters, one in each of the lower and middle redshift 
slices. For cluster B, the signs of deviation from the 
mean parameters in the core region in the lower redshift 
wall at z <~ 0.87 in contrast to the lack of deviation in 
the DS test of the central, overall supercluster redshift 
slice around z ~ 0.9 give further evidence that the cluster 
core lies in the higher z ~ 0.9 redshift slice since we ex- 
pect cluster core galaxies to show little deviation from the 
mean parameters. In the case of cluster C, the DS test on 
both redshift slices show signs of substructure emanating 
from the center, with the larger deviations present in the 
lower redshift slice. In 35.21 we motivated the choice of 
the z ~ 0.9 peak in the cluster C plot of Figure [6] as the 
one corresponding to the cluster core given its agreement 
with both the spectroscopically confirmed BCG and the 
distribution of the red-sequence galaxies significance con- 
tours. The deviations seen in the core of cluster C in this 
redshift slice could therefore be an indication that there 
is a line-of-sight group falling tow ards the cluster core 
(jEinasto et ai.ll2012t|Phiknev et al.lll996| ). The following 
search for individual groups will explore this possibility 
in more detail. 



6.2. 
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The distribution of possible substructure in our su- 
percluster field found with the DS test lead us to at- 
tempt to recover the substructure in a more tangible 
way. We used a modified frie nds-of-friends (FOF) algo- 
rithm (jHuchra fc Gellerlll982D with linking lengths of D L 
and Vl, representing distance and velocity separations, 
to identify associations of three or more linked galax- 
ies. Briefly, the algorithm uses each of our spectroscopi- 
cally confirmed cluster members and searches for friends 
within the two linking lengths. The test galaxy and its 
friends are removed from the main source list, and each 
friend is then tested against the remaining galaxies in the 
list. This is done iterative ly until there are no new friends 
found. If the final candidate friends list has three or more 
members then it is considered a FOF association, or else 
the galaxies are returned to the isolated galaxy list. We 
chose three member associations to compensate for the 
incompleteness of the spectroscopic data in an effort to 
identify regions that may be undersampled spectroscop- 
ically but still contain signs of structure. Monte Carlo 
Simulations are performed in ^6.31 to help determine the 
membership number above which we believe the FOF 
associations are tracing real structures. 

The linking lengths were determined through trial and 
error, with the aim of recovering structure while not link- 
ing all three cluster cores together. Monte Carlo Sim- 
ulations f £16.3[) were also used to test the FOF linking 
lengths, again to ensure that linking lengths were large 
enough to recover structure while not linking broad red- 
shifts walls across the test region. A more thorough 
method of determining the linking lengths would have 
been to construct a mock catalog from N-body simula- 
tions to test our linking lengths based on completen ess 
and recovery purity, as done by iFarrens et al.l ()2011[ ) in 
their search for groups and clusters in the 2SLAQ cat- 
alog (2dF-SDSS Large Red Galaxy and QSO survey). 
However, a search of the existing Millenium Simula tion 
light cones (|Lemson &: the Virgo Consortium! I2006D for 
a similar structure to the RCS 2319+00 supercluster did 
not yield any comparable structures and attempting to 
model this superstructure with our own simulations falls 
outside the scope of this paper. 

We found that the spatial linking length, Dl, was more 
influential in linking all three cluster cores as each clus- 
ter is separated by <3 Mpc from its nearest neighbor and 
Dl was thus limited to 0.5 Mpc. The velocity separa- 
tion, Vl, had more subtle effects on the group recovery 
rate and membership. We tested several different ve- 
locity separations, from 600-1500 km s _1 (~ 2-5x our 
spectroscopic redshift error of Az = 0.002). For the low- 
est velocity separations, the core groups of the clusters 
are recovered as several separate, small groups, even for 
members of the prominent peak in cluster C. We found 
that in our real data and Monte Carlo tests ( §6.3|l the dif- 
ference in group recovery between Vl lengths of 1000 and 
1500 km s _1 were minimal, recovering the same groups 
on average, with only one or two members missing from 
two of the larger groups. Since we suspect recent merger 
activity is present in some regions of the supercluster 
field and therefore expect our substructures may not be 
virialized, we choose to use the larger linking length of 
1500 km s _1 to look for extended structures. 

Figure [9] shows the results of the friends-of-friends test 
using the linking lengths of 0.5 Mpc and 1500 km s _1 , 



with linked FOF associations identified by same colored 
boxes. The red and blue lines on the plot indicate all 
supercluster member positions with the color and length 
of the lines representative of the relative velocity of the 
galaxies to a median redshift of z=0.897. The lines are 
oriented arbitrarily towards a center position on the plot, 
with the galaxy located at the end pointing towards the 
center for blue-shifted galaxies and at the end point- 
ing away from the center for red-shifted galaxies. This 
was done to retain some redshift information in order 
to better visualize the layout of the FOF associations 
within supercluster field. The grey contours represent 
the smoothed spectroscopic target density relative to 
the full optical catalog, with the higher target densities 
having darker, thicker contours. Due to an incomplete 
spectroscopic sample, some of our isolated superstruc- 
ture 'field' galaxies may in fact belong to infalling groups 
or filaments. We note, however, that FOF associations 
are still found in regions of lower spectroscopic density 
and some high target density regions contain few clus- 
ter members and no FOF associations. As such we do 
not suspect that our recovered structure is simply tracing 
the areas of high spectroscopic targeting. We identify 14 
FOF associations, with 7 having more than three mem- 
bers. Figure QJJ shows the number of groups for each 
membership number (hatched histogram). 

6.3. Monte Carlo Simulations 

To test the veracity of our group finding method and 
look for spurious groups in our identified FOF associa- 
tions we performed Monte Carlo simulations. All 1961 
confident galaxy rcdshifts were first randomly shuffled 
between their positions and then the new supercluster 
field members in the redshift range of 0.858 < z < 0.946 
were selected, now with their new positions. We then 
performed the same friends-of-friends test as done on our 
non-shuffled data. This was repeated 1000 times and 
the average number of FOF associations and member- 
ship number per association was compared with our real 
data. 

The average number of FOF associations returned in 
the MC tests was 8.8, compared to the 14 recovered over- 
densities in the real data. Figure [TUJ shows the number 
of FOF associations in each membership bin for the real 
data (red hatched histogram) versus the average numbers 
over all shuffled MC tests (black outlined histogram). 
The bins containing the FOF associations identified as 
the main cluster cores (see §6.41 for details) are labeled 
with their corresponding cluster ID. This shows that 
larger FOF associations are recovered only in the real 
data, including the cores of clusters C and A at 13 and 23 
members, but not in the randomly shuffled data, leading 
us to believe that FOF associations with > 5 members 
are tracing real group and core structures. Also, since 
the positions of the spectroscopic sampling are consis- 
tent between the real data and the MC test data, we 
are not recovering large structures simply due to highly 
targeted regions. 

6.4. Recovered Structures in the RCS 2319+00 
Supercluster Field 

Taking into account the results of the MC simulations, 
we designate all FOF associations with > 5 members re- 
covered from the FOF test on our real data as a lower 
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Figure 9. Friends-of-friends associations identified with the linking lengths of 0.5 Mpc and 1500 km s _1 are surrounded by 
like-colored boxes. Red and blue arrows represent red/blue shifted velocities of spectroscopically identified cluster members 
relative to a midpoint redshift of 0.897 in the 0.858 < z < 0.946 range and directed towards an arbitrary center position. Dotted 
circles represent 1 Mpc radii around the X-ray centres of the three cluster cores. Grey contours show the spectroscopic target 
density relative to the full optical source catalog, with higher target densities having thicker, darker grey contours (i.e. the ~0.5 
Mpc region around cluster A is a high target density region). 



estimate on the true substructure of the RCS 2319+00 
supercluster, noting that the simulations show the possi- 
bility that one of the recovered 5-8 member groups may 
still be a spurious detection (see Figure ITU|) . Several 
structural features are immediately discernible in Fig- 
ure [S] The cluster cores all comprise groups that seem to 
extend in an elongated fashion pointing to the neighbour- 
ing core, with 23 and 13 members in the core groups in 
clusters A and C. The galaxies in the core group of clus- 
ter A (Group 2) span a redshift range of ^6000 km s _1 
and have a dispersion of 1404 ± 249 km s _1 , with no in- 
terloper clipping. This is slightly higher, though within 
errors, than the dispersion we calculated in £15.31 from 
the spatial distribution and interloper clipping. In clus- 
ter C, the main core group (Group 10) contains galaxies 



spanning ^960 km s _1 with an estimated dispersion of 
722 ±182 km s _1 , very close to the dispersion calculated 
in ^5.31 with interloper clipping. 

Figure [TT] shows the relative velocity of the galaxies 
in groups from the average redshift of cluster A as a 
function of the spatial distance from cluster A, with the 
identified > 5 member groups highlighted against all con- 
firmed members. The numbers and colors plotted for 
each group correspond to the same numbers and colors 
as Figure [9j The top panel highlights the groups con- 
sistent with the spatial centers of the cluster cores and 
shows that the core group of cluster A is elongated in ve- 
locity space as well as in its spatial distribution, while the 
core group of cluster C is more tightly bound. The lower 
redshift galaxies in the core group of cluster A may be 
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Figure 10. The number of FOF associations found for each 
membership bin for the real supercluster data (red hatched 
histogram) and the MC test data described in %.3\ (aver- 
aged over all tests - black outlined histogram) recovered using 
the FOF method. The letters indicate the bins in which the 
groups corresponding to the three cluster cores are found. 

interlopers from the lower redshift structure, from which 
they seem to extend. 

Once again, cluster B presents some difficulties. The 
two redshift peaks that were apparent in Figure[5]fall into 
two separate groups with 8 members at an average red- 
shift of 0.870 and 4 members in a group ^5700 km s _1 
away at z=0.906. Again, we suspect that more galax- 
ies belong to the higher redshift group in this core, cor- 
responding to the red-sequence density, X-ray contours 
and BCG, that have yet to be targeted with spectroscopy. 
Figure [TT] again shows that the smaller group in cluster 
B (group 9) lies along the same velocity space as clusters 
A and B. It can also be seen in Figure [9] that this higher 
redshift group falls in the region of lower spectroscopic 
target density. In light of this, we assign the 4 member 
group 9 as the core group of cluster B and the 8 member 
group 6 is considered a potentially infalling group dis- 
tinct from the core of the cluster. We attempt to correct 
for spectroscopic completeness in Section [51 

An interesting structure recovered through the FOF 
algorithm is an apparent 9 member filamentary struc- 
ture (Group 5) stretching from cluster B towards cluster 
A, first hinted at in the results of the DS test around 
the z ~ 0.90 redshift wall (S [671] Figure |8j. This struc- 
ture has an average redshift of z = 0.906 with galaxies 
spanning ^2800 km s _1 towards higher redshift (bot- 
tom panel, Figure [TT1) . It corresponds well to the pre- 
viously confirmed H erschel infrared bright filament that 
iCoppin et al.l (|2012f) found may be forming stars with a 
SFR ~ 900 M Q yr _1 and building up significant stellar 
mass. 

Other notable groups include a 6 member group cen- 
tered between the three component cluster cores (group 
7) as well as an 8 member filamentary group falling in to- 
wards cluster C (group 13). Both of these groups appear 
to be falling into the clusters from the lower redshift wall 
(bottom panel, Figure ITTj) and may represent infalling 
groups. 

Comparing our recovered groups from the FOF method 



with the substructure detected with the DS test (see the 
DS test bubble diagram, Figure [8]), we see that the cen- 
tral group (7), the filamentary structure (group 5), and 
the 8 member group in the line of sight of cluster B 
(group 6) correspond to regions of deviation in the DS 
test. Cluster C has one large core group (10), which is 
relatively compact in velocity space, and a small, three 
member FOF association (12) found in the same line- 
of-sight. This small overdensity falls within the lower 
redshift wall and accounts for the deviation seen in the 
DS test of this redshift slice when tied with the other 
blue-shifted galaxies seen south-west of the cluster core 
in Figure [SJ This area and other areas of possible sub- 
structure in the DS test found to have small, 3-4 mem- 
ber overdensities with the FOF method (groups 4 and 8), 
lack adequate members in the immediate area to confirm 
substructures. The addition of further spectroscopic red- 
shifts may serve to confirm the DS substructures. 

We find that of our 302 cluster members in the super- 
cluster field ~13% are found in the core cluster groups 
and ~18% belong to the potential infalling groups. We 
note that these are rough estimates, since our superclus- 
ter 'field' galaxies may in fact be tied into groups given 
further spectroscopy. Areas of deviation from the global 
properties in the DS test that are not recovered with 
the FOF method may belong to structures whose linking 
members have not yet been confirmed spectroscopically, 
leaving them out of the FOF group list and identified as 
part of our superstructure 'field' population. Conversely, 
we may find a higher proportion of iso lated galaxie s given 
further spectroscopy. At low redshift ISmall et all (1998) 
found that approximately one third of the galaxies associ- 
ated with the Corona Borealis supercluster do not belong 
to the component clusters. At higher redshifts we would 
expect the number of galaxies in superclusters without 
specific cluster associations to increase as the c lusters are 
less re laxed. This has already been seen by iGal et all 
(2008) in their study of the structure of the CL1604 su- 
percluster at z ~ 0.9, where they found that a large 
portion of galaxies were associated with the connecting 
large-scale filaments and not with the individual clus- 
ters and groups in their supercluster. As such, while we 
do expect that additional spectroscopically confirmed su- 
percluster members will result in more discernible struc- 
tures, the percentage of unassociated supercluster mem- 
bers will most likely remain high. 

7. THE COLOR-MAGNITUDE RELATION 

A color-magnitude relation (CMR ) for e lliptical field 
galaxies was first identified by iBauml (|1959D , with bright 
ellipticals having redder colors then their fainter coun- 
terparts. In clusters, the CMR for early-type galaxies, 
the red-sequence, has been sho wn to have a tight, well- 
defin e d slope at low redsh i ft (eelVisvanathan fe Sandaeel 
[19771 iBower et alj U99l lLopez-Cruz et al.l 120041) and 
the cluster cor es are dominated by early-type galaxies 
(Dressier 1980). At higher redshifts however, the dynam- 
ically younger clusters have shown an increasing scatter 
in their CMR and a deficit of red-sequ ence galaxies at 



the lo w- luminosity end o f the slope (e.g. [De Lucia et al 
20041: iTanaka et at 2005; Mei et alj[2"009t iLemaux et al 



20121) . This is seen as a sign that the cluster red- 



seque nce is still in t he pr ocess of being built up at 
z~l. iRudnick et all ([2012) propose that the lack of low- 
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Figure 11. Velocity versus radius for all confirmed cluster 
members in the RCS 2319+00 supercluster relative to cluster 
A. Top panel: Groups that correspond to the spatial positions 
of the three cluster cores are highlighted: cluster A - group 
2, cluster B - groups 6 and 9, cluster C - group fO. The two 
groups near the spatial center of cluster B lie on two of the 
distinct redshift walls that extend through the supercluster 
region. We assign group 4 as the cluster core as it falls along 
the same velocity space as the two other cluster cores and 
corresponds to the red-sequence peak redshift in this region. 
Bottom panel: Groups having >5 members outside of the 
cluster cores are highlighted. Two larger groups are seen to 
be potentially falling into the cluster from the lower redshift 
peak. The filamentary structure between clusters B and C 
(group 5) is shown to lie within the same redshift wall as the 
core groups. 

luminosity red-sequence galaxies may be due to the ac- 
cretion of faint blue galaxies whose star-formation is sup- 
pressed in groups or in the cluster and experience mergers 
with other red-sequence galaxies, increasing their mass 
and moving them along the red-sequence. Recent stud- 
ies of clusters a z > 1 differ over the basis for the scatter 
in the red-sequence at higher redshifts, with clusters at 
z = 1.46 and 1.62 found to be cons istent with models 
of a single burst of star- formation (Hi lton et al.1 120091 : 
iPapovich et al.l 120101 respectively) while more diverse 
star formation histories are needed to expl ain the width 
of th e CMR in clusters at 1.0 < z < 1.5 (jSnvder et all 
I20T1 . 

7.1. Color-magnitude diagram 

In this section we investigate the color-magnitude rela- 
tion in the RCS 2319+00 supercluster. The R c -z' versus 
z' color-magnitude diagram for all spectroscopically con- 
firmed galaxies within the supercluster redshift range of 
0.858 < z < 0.946 is shown in Fi gure [T2l with the red- 
sequence slope at z = 0.9 (from IGilbank et al.l [2008b) 
plotted for reference as a solid line, and dotted lines de- 
noting ±0.2 mag from the the red-sequence. The errors 
on the photometric catalog range from 0.01 mag for the 
brightest galaxies, increasing to 0.62 mag for the faintest, 
z 1 = 24.85 magnitude cluster member and 0.73 mag for 
the faintest R c = 25.73 magnitude member, with an av- 
erage error of ~ 0.08 mag in both filters. 

The two top panels of Figure [T3] both show all spectro- 
scopically confirmed RCS 2319+00 supercluster members 
plotted as black circles. The colored symbols in each 
panel represent members of the identified groups from 
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Figure 12. R c — z' versus z color-magnitude diagram. Black 
points represent all spectroscopic galaxies in the redshift 
range 0.858 < z < 0.946. The top two panels plot all members 
of the RCS 2319+00 structure with colored symbols highlight- 
ing galaxies belonging to the numbered groups identified with 
the FOF algorithm as shown in Figure[9] divided into the core 
groups (top panel) and > 5 member groups (middle panel). 
The bottom panel shows the comparison field sample from 
the large RCS IMACS spectroscopic campaign (R. Yan et 
al. 2012, in preparation). Lines represent the red-seq uence 
slope at z=0.9 (solid line, from IGilbank et al.l I2008bf) and 
±0.2 of the red-sequence slope (dotted line). Note that no 
correction for redshift success has been applied to the color- 
magntude diagram. 

our friends-of-friends algorithm in with the numbers 
and colors plotted for each group corresponding to the 
same numbers and colors on the structure plot (Figure[S]) 
separated into the two panels to highlight the cluster core 
groups (top panel) and the > 5 member groups (middle 
panel) . The bottom panel shows the color- magnitude di- 
agram of a comparison field sample, detailed below. It 
should be noted that the redshift success rate of the spec- 
troscopic sample is a function of both the magnitude and 
color for galaxies fainter than z' ~ 22 mag (see FigureS]), 
and is not accounted for in the color-magnitude diagram 

Due to the incompleteness of our spectroscopic sample, 
we attempt a rough completeness correction when look- 
ing at the overall population distribution in the color- 
magnitude diagram. Since the spectroscopic correction 
weights do not take into account the color of the galaxies, 
we instead use the redshift success rates for bright (z' < 
22 mag) and faint (z' > 22 mag) galaxies found in the 
right panel of Figure [4] We set the success rate of the 
bright galaxies conservatively at 70% from our 70-80% 
redshift success rate for bright galaxies. For the faint 
galaxies, we define a color cut at R c — z' = 1 and find 
that bluer galaxies have a success of ^50% while the red- 
der galaxies have a redshift success of only ^25%. We 
can use these rates to correct the faint galaxy numbers 
in Figure [12] up to the redshift success rate of the bright 
galaxies in order to better compare the populations. This 
is an oversimplification of the completeness correction, 
however it is sufficient given our small number statistics. 

Using all spectroscopic structure members corrected 
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to the redshift success of bright galaxies and defining a 
rough limit for the red versus blue galaxies at -0.2 mag 
of the red-sequence line, we find 36±3% of the galaxies 
are red and 64±4% of the structure members are blue 
galaxies, with Poisson errors. We can compare these 
populations to the field in the same redshift range us- 
ing a spectroscopic sample of galaxies derived from ob- 
servations whose target cluster falls outside this redshift 
range in the RCS IMACS spectroscopic campaign (R. 
Yan et al. 2012, in preparation). We find 422 spectro- 
scopically confirmed field galaxies that fit this criteria 
in the IMACS sample. The field sample is plotted in 
the bottom panel of Figure [T^] with the red-sequence line 
again at z — 0.9 again plotted for reference. Those galax- 
ies populating the red-sequence make up only 18±2% of 
the sample while the blue cloud is populated by 82±4% 
of the comparison field sample. The limit on the depth 
of the spectroscopy is brighter for the field sample at 
z' <~ 24 mag (M z > ~ —18.35) compared to the limit of the 

RCS 2319+00 sample at z' ~ 25 mag (M z > 17.35). 

As such, we can set the depth of the RCS 2319+00 struc- 
ture to that of the field sample for better comparison. 
The fraction of red and blue galaxies in the structure to 
this brighter depth remains unchanged. One caveat to 
the field sample is that, though it is drawn from spectro- 
scopic observations whose main target cluster lies out- 
side of the supercluster range, there may exist groups 
and low-richness clusters within the supercluster range 
that have not yet been identified. Despite this, our su- 
percluster galaxies show an excess red-sequence over the 
field population, as expected in clusters. To further in- 
vestigate this excess we look at the identified structure 
from $6l 

It can be seen that the cluster cores (groups 2, 9, 
and 10) show a more prominent red-sequence than the 
groups, with the core groups having 77+13% red galax- 
ies and only 23+7% in the blue cloud, for all identified 
members down to the fainter z' ~ 25 mag limit. The 
large errors are due to the small numbers, even after 
correcting for the redshift success. The groups show a 
lower fraction of red-sequence galaxies at 37+10%, with 
the percentage of blue cloud galaxies almost three times 
higher than those in the cores at 63+13%. Though we 
do not expect to see significant differences between the 
populations given the small numbers of galaxies in our 
sample and the large errors, the populations are still not 
consistent within errors, only crossing at ~ 2a. From this 
we infer that the two group categories contain galaxy 
populations that may be in slightly different stages of 
evolutio n. In their study of the C11604 supercluster at 
z ~ 0.9, iLemaux et al.l (|2012f) found that the fraction of 
their cluster and group galaxies in the red-sequence was 
identical at 47%, suggesting a similar level of processing 
in the two environments. They also found that 70% of 
the red-sequence galaxies in their supercluster belong to 
the clusters and groups. Combined, the core and group 
members of the RCS 2319+00 supercluster account for 
only 38+5% of the corrected red-sequence supercluster 
members, which again hints that some fraction of the iso- 
lated supercluster 'field' galaxies may belong to larger, as 
yet unidentified structures. In the supercluster isolated 
'field' population, which accounts for 226 of our 302 su- 
percluster members, 29+3% of the corrected galaxies are 
red while the remaining 71+5% are in the blue cloud 



region, consistent with the group fractions. 

The interesting filamentary structure between clusters 
A and B (group 5) has 8 of its 9 identified members in 
the blue cloud and includes some of the faintest identi- 
fied members so far, two of whom are also among the 
bluest galaxies currently identi fied in the overall struc- 
ture. This is not surprising as iCoppin et all (|2012[ ) have 
shown that this region is undergoing an excess of star for- 
mation activity. If we remove this unique structure from 
our group category, the remaining >5 member groups 
contain 43+12% red-sequence galaxies, still lower than 
the fraction of red-sequence galaxies in the cluster cores, 
though now consistent within < 1.5<7. This slightly 
higher gr oup red-fraction is in closer agreement with the 
results of lLemaux et al.l (|2012l ) for both their groups and 
cluster galaxies. On average then, we believe the groups 
in our structure have undergone some processing, though 
less than the cluster core galaxies, with the filamentary 
structure standing out as a unique environment of blue, 
star-forming galaxies currently undergoing intense activ- 
ity. 

Recalling that the Monte Carlo simulations run in ^ 36.31 
showed the possibility of one spurious 5-8 member FOF 
group, we look to sec what effect this might have on the 
red-sequence fractions of the populations. The validity 
of the core groups is supported by the optical color data 
and the X-ray detections and we are therefore confident 
that the FOF overdensities represent real groups. The fil- 
amentary structure is detected in the infrared and again 
we assume that it's corresponding FOF group can there- 
fore be assumed to trace a real structure. This leaves 
only one of the four remaining > 5 member groups as 
the possible spurious structure detection. If we remove 
the galaxies belonging to a single group and reassess the 
fraction of red vs blue galaxies remaining for the group 
population, repeating the process until each group has 
been removed once, we find that the fraction of red and 
blue galaxies does not change significantly, with the frac- 
tion of red galaxies varying from 41-45% with errors of 
< 14%, consistent with the previous results. 

With the aid of the color information, we can revisit 
the assignment of group 9 as the core group of cluster 
B over group 6, the lower redshift group also consistent 
with the spatial position of the red-sequence overdensity 
and the X-ray emission. In Figure [T2l it can be seen that 
all of the galaxies currently identified in group 9 are red, 
while 50% of the identified galaxies in group 6 are blue, 
closer to the fraction seen for the groups. This lends more 
credence to our claim that group 9 is the core group of 
cluster B while group 6 is a possible infalling group from 
the lower redshift wall. 

7.2. Color Distribution In Different Environments 

To further investigate the populations in the different 
environments and look for possible signs of the build- 
up of the red-sequence we look to the distribution of 
galaxy colors in the different density regions. Here we 
do not attempt to correct for the redshift success based 
on magnitude or for the galaxies mass or luminosity dif- 
ferences, which will be the subject of a follow up paper. 
Instead, we normalize the number counts in the differ- 
ent populations to one in order to look at the relative 
color distributions. In Figure [131 we plot the normalized 
number counts of galaxies in the different environments, 
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Figure 13. Number of galaxies in each environment normal- 
ized to one versus R c — z' color for the spectroscopically con- 
firmed supercluster members plotted as a difference over the 
normalized comparison field sample. Errors given are Pois- 
son errors found for each population. The solid line shows the 
color distribution in the cluster core galaxies over the field, 
the dashed line represents the > 5 member, non-core groups 
with the filament removed, the dash-dot line shows the fila- 
ment galaxies and the dash-dot-dot-dot line encompasses all 
remaining isolated supercluster 'field' members. 

binned in color and plotted as a difference over the nor- 
malized comparison field sample discussed above. We 
first k-correct the apparent magnitudes using the spec- 
troscopic redshifts and ve rsion 4.2 of the kcorrect code by 
iBlanton fc Roweisl (|2007| ) to better compare the samples. 
The solid line shows the color distribution in the cluster 
core galaxies over the field, the dashed line represents 
the > 5 member, non-core groups with the filament re- 
moved, the dash-dot line shows the filament galaxies and 
the dash- dot- dot- dot line encompasses all remaining iso- 
lated supercluster 'field' members. The Poisson errors 
plotted are based on the parent population only, with no 
scaling applied when the difference between each popu- 
lation with the field population is found. 

The isolated supercluster galaxies possess the expected 
bimodal distribution of blue and red galaxies over the 
field, with a distinct, underpopulated 'green valley', 
though ~ 38% of our spectroscopic cluster galaxies fall in 
this intermediate color area. The core and group galax- 
ies both show peaks in the red end of the histogram, as 
expected, with few blue galaxies and a deficit of 'green' 
galaxies compared to the field. The group galaxies do 
show a slight excess of galaxies at the bluest end of the 
distribution, though the large errors render the peak con- 
sistent with the comparison field at this color. 

The filament between clusters A and B shows three 
peaks above the field: a small peak at the red sequence 
end; one populating the green valley; and a final excess at 
the bluest end of the color distribution. The small num- 
ber of filament galaxies currently in our spectroscopic 
sample (9 members) results in large errors and there- 
fore the three filament peaks are consistent with the field 
within the errors. Though we do not have statistically 
significant data with such small sample numbers in the 
green galaxy population from the filament, we propose 
that, if this peak is real, these intermediate color galax- 



ies may in fact belong to a transient population in the 
process of migrating to the r ed-sequence, similar to those 
seen bv lBalogh et al.l ([2011ft in their study of seven spec- 
troscopically confirmed galaxy groups at 0.85 < z < 1. 
They found that 30% of their confirmed group members 
populated this intermediate 'green' color region and pro- 
pose that these galaxies represent a transient popula- 
tion who have recently ceased their star formation and 
are in the midst of migrating from the blue cloud to 
the red-sequence. Currently, two thirds of the filament 
galaxies are found in this intermediate area of the color- 
magnitude diagram versus the ~ 45% of field galaxies 
at this color. Further studies with multiwavelength data 
on the RCS 2319+00 supercluster is needed to determine 
the properties of these galaxies. 

8. LUMINOSITY FUNCTIONS 

We wish to put a rough estimate on the average mass 
of our identified groups in order to better trace the 
buildup of our supercluster and to place our groups in 
context with theories o f hierarchical f ormat i on from semi- 
analytic models (eg. IMcGee et all 120091 : iBerrier et al.l 
I2009t iDe Lucia et all 1201% iCohd 12012ft . The majority 
of our groups have too few spectroscopic members to 
compute reliable velocity dispersions and, as we have 
seen in M5.41 masses derived from velocity dispersions can 
carry large errors even with a larger sample. Instead, we 
attempt to estimate the mean mass of our groups by 
fitting luminosity functions (LF) to two stacked group 
categories; the cluster cores and the potential inf ailing 
groups. 

Luminosity functions measure the number of galaxies 
per volume in a given luminosity bin and can be used 
to test theories on the formation of galaxies, their evo- 
lution, and what role the environment plays on these 
factors. The faint-end of the luminosity function is pop- 
ulated by a mix of faint red and blue galaxies and slight 
variations in its slope can be used to help distinguish be- 
tween different structure formation theories. A steeper 
faint-end slope may be indicative of the accretion of faint 
galaxies onto the group or cluster while a flatter slope 
would be expected in a relaxed cluster dominated by 
early - type galaxies (eg . IGoto et al.l 120021: iPropris et al.l 
[20031: IGoto et al.l [20051: IZenteno et all 12011ft . LFs can 
also be used to estimate the number of galaxies in each 
group or cluster by integrating the function to given lu- 
minosity limits. While cluster luminosity functions are 
generally constructed from background subtracted pho- 
tometric catalogs, we attempt to measure the luminosity 
functions of the stacked clusters and groups in our super- 
cluster through the recovered spectroscopic membership 
from the friends-of-friends search for substructure in the 
previous section. 

We use the fu nctional form of the Schechter function 
(Schcchterl ll976ft to fit the luminosity function, 

<MM) = 0.41n(10)^(10°' 4(M *- M) ) 1+Q exp(-10°' 4(M *- M) ), 

(6) 

where a is the slope of the faint-end power-law, M* is 
the characteristic magnitude and <fi* is the normalization 
of the Schechter function. 

8.1. Luminosity Function of the RCS 2319+00 
Supercluster Field 
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Figure 14. LF of the best-fit Schechter function to all com- 
pleteness corrected spectroscopic supercluster members in the 
35 x 35 arcminute region with the BCGs from the three clus- 
ter cores removed so as not to influence the bright-end fit of 
the LF. 

We first construct a LF for the full catalog of 302 
spectroscopically confirmed supercluster members in our 
35 x 35 arcminute region to estimate the value of M*, and 
look at the shape of the faint-end slope. First, we con- 
vert the photometry of all cluster members to absolute 
magnitude M z > using, 
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Dl(z) 



10 pc 



K(z) 



(7) 



where Dl(z) is the luminosity distance to our aver- 
age supercluster redshift and K(z) is the k-correction 
computed using version 4.2 of the kcorrect code by 
iBlanton fc Roweii (|2007| ). Our apparent magnitudes 
(m z <) have already been corrected for galactic dust ex- 
tinction. Th e spectroscopic magnitude weights calcu- 
lated in ^4.21 are applied to each confirmed member with 
no geometric correction as we are looking at the statistics 
over the whole field and not at individual galaxy prop- 
erties. The completeness weighted spectroscopic mem- 
bers were then binned by absolute magnitude in 0.25 
mag bins and the bins for each core or > 5 member 
group were stacked. Fractional errors were added to 
the stacked, weighted number counts in each magnitude 
bin. The Schechter function was fit using the Levenberg- 
Marquardt algorithm for non-li near least squares fitting 
using the mpfit routine for IDL ()Markwardtll2009D . 

Figure Q3] shows the LF for all weighted, spectroscopi- 
cally confirmed members in the RCS 2319+00 superclus- 
ter. The BCGs from each of the three cluster cores 
were removed from the fitting so as not to influence the 
bright-end fit of the LF. The solid line shows the best fit 
Schechter function. The characteristic magnitude M*, is 
estimated to be M*, = -21.01 ± 0.22. This value is 
slightly lower than the value derived by Gilbank et al. 
(2008bJ) of M*, = -21.39 at z = 0.9 using the passively 



evolving model on a fit to the LF of composite RCS clus- 
ters at z = 0.4. We note that our fitted value of the faint- 
end slope of a = —1.30 ± 0.17 is at the high-end of the 
values found with other galaxy clusters in the literature. 
The fit to the composite RCS clusters of iGilbank et al.1 
(20083) re sulted in a faint-end slope of a = -0.94±0.04. 
iLin et all (|2004[ ) used a composite K-band LF of 93 X- 
ray selected groups and clusters from the Two Micron All 
Sky Survey and found a faint-end slope is fit by a slope 
ranging as —0.84 < a < —1.1. The closest agreemen t 
to our estimate of a was found by iZenteno et al.l (|201lD , 
who measured an average faint-end slope of a rs —1.2 
for their four Sunyaev-Zel'Dovich selected clusters in the 
r, g, i, and z-bands. 

To construct our supercluster LF, we use the 
weighted number counts of all spectroscopically con- 
firmed cluster members over the broad supercluster re- 
gion rather than using background subtracted number 
counts from pho tometry in the much smaller rad ii of r^oo 
Lin et al.ll200l. 0.5 r9nn (|Gilbank et al.ll2008bD or r 200 



Zenteno et al.ll201ll ) around groups and cluster cores as 



done in the previous studies. The steeper than average 
faint-end slope we find for RCS 2319+00 may therefore 
be reflective of the presence of infalling structures and 
individual galaxies in the lower density regions of the 
supercluster field. The addition of these lower density 
regions to our LF may also be the cause of the slightly 
lower charact eristic magnitude foun d here compared to 
that found bv IGilbank et all (|2008bQ . However, we must 
again note that our fit is based on an incomplete, com- 
pleteness corrected spectroscopic sample, which may also 
affect the shape of the LF. For these reasons, we choose 
to use the values for the M*, and a par ameters found 
for the composite RCS galaxy clusters by IGilbank et al.l 
(2008b]) when looking at the LFs of the core and > 5 
member groups in the following section as they appear 
to be more representative of the denser cluster cores and 
groups than our fitted values from the broad supercluster 
field. 

8.2. Luminosity Functions of the Cluster Cores and 

Groups 

To construct the LFs of our stacked > 5 member 
groups and the cluster cores recovered from the friends- 
of-friends output we follow the same methods described 
above for the whole supercluster member catalog with 
the exception of the spectroscopic completeness weights, 
which were calculated using the photometric catalog in 
a smaller radius around the mean group position (see 
^4. 21 for details on the spectroscopic weight calculations) 
rather than the completeness weight for the full 35x35 
arcminute field of view. This was done to simulate a 
smaller group radius of ~ r2oo , which we initially set at 
1.5 Mpc for the cluster cor es, based on the valu es found 
for the clusters using the lEvrard et alj (|2008l ) relation 
(see Table [5]) and an arbitrary radius of 1 Mpc for the 
> 5 member groups. Since the LFs of our recovered 
groups and cluster cores are based on a small number of 
completeness weighted spectroscopic sources per group, 
we find that fitting all three Schechter parameters simul- 
taneously is not reliable. Instead, we fix the characteris- 
tic magnitude and faint-end sl ope to those measured for 
the composite RCS clusters bv IGilbank et al.l (|2008bD of 
Ml, = -21.39 and a = -0.94, as motivated above. For 
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the stacked LF of the three core groups (2, 9 and 10) the 
BCGs were again removed from the fitting so as not to 
influence the bright-end fit. 

The stacked luminosity functions are shown in Fig- 
ure [15l The LF of the stacked cluster cores is shown 
as closed circles with the best-fit function plotted as a 
solid line. The open circles and dotted line show the re- 
sults for the stacked > 5 member groups. The dashed 
vertical line indicates the position of Al*, for reference. 
The luminosity functions with fixed Al*, and a were fit to 
<f>* , integrated to estimate the number of galaxies in the 
stacked groups, N ga i s , and then divided by the number 
of groups in the stacked LF as a proxy for the mean rich- 
ness of the groups. We used two methods of estimating 
the group mass from the number of galaxies. 

The first method employs the halo occupation num- 
ber (HON), based on the halo occupation distribution 
(HOD); the probability distrib ution of finding N galax- 

BensonTl200ll: 
2006!). In the 



ies in a halo of mass Al (eg. [S eljak 



iBerlind fc Weinberg 2002; iLin et al 



2000; 



2004, 



HOD formalism, the HON is defined as (N)(M), the 
mean number of galaxie s per halo , and c an be estimated 
by integrating the LF. ILin et aLl (|2004l ) used a sample 
of 93 X-ray selected groups and clusters in the nearby 
universe to form a composite K-band LF and found a 
N 500 — Mpon scalin g relat ion. Extending their work to 
z = 0.9, ILin et all (I2006T) found no signs of evolution 
in the scaling relation. iZenteno et al.l (|2011h converted 
this relation to -/V200 — M200 f° r their work with Sunyaev- 
ZeTDovich selected clusters and found that the ILin et akl 
(2004) relation becomes: 



N> 



200 



(36 ± 3) 



M: 



200 



10 14 /ifo 1 M G 



0.87±0.04 



(8) 



The N — M relation used is based on the number of 
galaxies within the viral radius of the group or cluster. 
Since we do not have that information available and our 
LFs are constructed from weighted spectroscopic samples 
with small numbers of known galaxies, we retain our 1.5 
and 1 Mpc completeness corrections initially for these 
estimates for the cores and > 5 member groups respec- 
tively. The total number of galaxies A/200 is estimated 
from the integration of the luminosity function. Here 
the integration is done to a limit of M* +2. 

Using the estimated A/200 and its resulting mass, a new 
estimate of the virial radius r2oo can be found and a new 
spectroscopic weighting for the galaxies can be performed 
to reflect this new radius. This is done iteratively until 
the resulting radius matches the input radius. The mass 
estimates with the initial and final radii vary by < 0.1 x 
10 14 M©, with the core mass decreasing slightly and the 
group mass increasing for the final, smaller radii. For 
the cluster cores, we find an average N ga [ s = 23.98 and 
a mas s estimate of ~ 0.9 x 10 14 M Q using the ILin et al.l 
(2004) A/200 — M200 scaling relation with the galaxies 
completeness weighted to an average virial radius of ~ 
0.7 Mpc. The mass estimated for the >5 member groups 
is over a factor of two lower at ~ 0.4 x 10 14 M© for N gals = 
11.55 per group in a ~ 0.5 Mpc radius. 

The second method we use to estimate the masses of 
our groups from the stacked LFs involves the cluster and 
group richness estimate used with the Sloan Digital Sky 
Survey maxBCG catalog (e.g. lHansen et al.ll2005l 120091 : 
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Figure 15. Stacked LFs for the three cluster cores in the su- 
percluster from 1.5 Mpc completeness corrected spectroscopic 
members found using the FOF algorithm (closed circles, solid 
line) and for the 1 Mpc completeness weighted members in the 
five > 5 member recovered groups (open circles, dotted line) 
in magnitude bins of ±0.5 mag. The curves show the best- 
fit Schechter functions us ing M*, = —21 . 39 (ve rtical dashed 
line) and a = -0.94 from lGilbank et all (|2008bf ). 



iBecker et all 120071; Uohnston et al.l 12007). In this case 
the richness, N200, represents the number of red galax- 
ies that are fainter than the BCG in an area that ap- 
proximates r2oo- N200 is first calculated for an arbitrary 
radius, which we take to be the 1.5 and 1 Mpc radii 
used above, and then calculated using the new r2oo from 
the maxBCG -/V200 ~ ^200 relations. Richness-size rela- 
tions have been found for the maxBC G catalog from both 
the radial density profiles of clusters (|Hansen et al.ll2005[ ) 
and weak lensing (jJohnston et al with the weak 

lensing relation scaling to the cluster halo mass. 

Due to our small numbers and completeness weighting, 
which is not dependant on color information, we can not 
readily extract the red and blue galaxies into separate 
LFs to find different N ga i s for the two populations. In- 
stead, we use the rough percentage of red-sequence galax- 
ies in each of the group categories found in ^ IZ. II of ^77% 
for the core groups and ~33% for the >5 member groups. 
We used these percentages to scale the N ga i s recovered 
from the LFs to an approximate A/200 to be used with the 
maxBCG relations. For the cores, we find A/200 ~ 41.2 
from r2oo ~ 1.1 Mpc and subsequent mass estimates of 
~6x 1O 14 M per cluster. For the >5 member groups 
we find r 2 oo ~ 0.6 Mpc and A/200 ~ 6.7, resulting in mass 
estimates of ~ 5 x 10 13 M Q per group. 

The average group mass estimates resulting from the 
different N—M scaling relations are relatively consistent, 
with the cluster cores having masses ~ 1O 14 M0, an order 
of magnitude more massive than the > 10 13 M Q estimates 
for the potential infalling groups. The mass estimate for 
the cluster cores of ~ 0.9 — 6 x 10 14 M Q is consistent with 
those calculated from the spectroscopic velocity disper- 
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sions and the X-ray data of 10 14 5 — 10 14 ' 9 M©. From this 
we reason that the average group mass from our LFs 
of ~ 4 — 5 x 10 13 M© is a reasonable initial group mass 
estimate. 

As a test, the LFs were also fit with only the faint- 
end slope set and allowing M*, to vary. This was done 
in order to check that any mass segregation between 
the denser cluster cores and the surrounding groups in 
the supercluster field did not cause the characteristic 
magnitude to shift drastically in the LFs of the stacked 
core/group galaxies (Figure [TJ] We found only slight 
changes in the value of M*, that did not result in any 
changes in the resulting average cluster core or group 
mass estimates. 

8.3. Comparison to Theories of Hierarchical Structure 
Formation 

In order to put our ~ 10 13 M© recovered groups in con- 
text, we look to studies which employ semi-analytic tech- 
niques to retrieve galaxy merger tree information from 
N-body simulations in order to trac e the merger history 
of cl usters. In two recent studies bv lBerrier et al.l (2009) 
and iMcGee et all (|2009f ). slightly differing results were 
obtained as to the fraction of galaxies in clusters today 
that were accreted onto the clus t er as a member of a 
~ 10 13 M group. IBerrier et all (|2009D found that, of 
galaxies residing today in 10 142 Mq clusters, < 17% 
were accreted in groups of > 5 members, and ~ 25% 
from halos wit h masses of ~ 10 13 fe " 1 Mq. For the same 
mass clusters, IMcGee et al.l (|2009l ) found a higher frac- 
tion of galaxies fell into the cluster from ~ 10 13 h~ Y Mq 
groups, at ~32%. The discrepancy in the results of the 
two studies is likely due to the different methods used to 
p opulate the d ark matter halos with galaxies. 

iCohnl (|2012f ) also found that ~ 30% of galaxies in clus- 
ters with mass > 10 14 h~ l Mq at z = 0.1 had been 
accreted onto the final cluster as members of halos with 
> 10 13 /i _1 M©, with the most massive cluster galaxies 
having a larger fraction coming from these groups. This 
work also found that the largest subgroups accreted onto 
clusters were aligned along the clusters major axis, which 
we see in Figure [§] in groups 5 (the filamentary group 
between clusters A and B), 6 (large group in the line-of- 
sight of cluster B) and 13 (the large filamentary group 
seemingly falling into cluster C). 

Another study, by IDe Lucia et al.l (|2012l ). also used 
galaxy merger trees to study the accretion of galaxies 
onto clusters. While their study focuses primarily on 
the stellar mass of the galaxies accreted through differ- 
ent mass groups, they also find that a large fraction of 
the most massive cluster galaxies and about half of the 
low to intermediate stellar mass cluster members fell into 
the final cluster from ~ 10 13 M© groups. These studies 
all suggest that some pre-processing of cluster galaxies 
in groups is occurring prior to accretion onto the final 
cluster. 

The three clusters in RSC2319+00 are thought to be 
on a collision cou rse to form a massive , > 10 15 M© clus- 
ter by z = 0.5 (jGilbank et all l2008ah . If we consider 
cluster A to be the most massive halo in our superclus- 
ter (based on the X-ray masses) and assign it as the final 
host halo in the eventual merging of the RCS 2319+00 
supercluster, the spectroscopic members (as opposed to 



the completeness corrected members) of the remaining 
~ 10 13 Mq groups and clusters currently account for 
22% of the spectroscopic cluster members. We recall 
that our other identified groups (< 5 members) as well 
as some of our isolated supercluster field galaxies may 
in fact be identified as larger groups given further spec- 
troscopy and, on the other hand, we may find a larger 
isolated galaxy population in the supercluster field so this 
percentage may change. 

While comparing to these studies gives a general idea 
that our recovered groups likely represent real substruc- 
ture in our supercluster that is destined to be accreted 
into the final cluster, a comparison to clusters at sim- 
ilar redshift a nd of similar final mass is more prudent. 
IMcGee et al.l fj2009f ) extended their study to trace the 
history of accretion onto galaxy groups and clusters at 
different redshifts. They determined that at z = 1.0 clus- 
ters with masses of ~ 10 14 5 h~ x Mq had ^36% of their 
galaxies accumulated through > 10 13 h~ x Mq halos. For 
Coma-like clusters of ~ 10 15 M© at z = 0.5, which 
the RCS 2319+00 supercluster is predicted to form, they 
found that ~ 40% of the cluster galaxies had been ac- 
creted through large halos over the merger history. Of 
more relevance is their result that by tracing the merger 
trees of the most massive M c ; us ~ 10 151 h^ 1 Mq galaxy 
clusters today, they were found to never have more than 
17% of their galaxies in halos of mass 10 13 — 10 14 Mq, 
which occurred at a lookback time of ~7 Gyrs, about 0.3 
Gyrs more recently than the epoch at which we are view- 
ing the RCS 2319+00 supercluster. In this context, our 
22% of galaxies in the supercluster currently residing in 
groups with masses of 10 1 3 — 10 14 M^, thoug h slightly 
high, is comparable to the Mc Gee et al.l (2009) findings. 
Again we reiterate our numbers are estimates based on 
incomplete spectroscopic coverage. 

9. SUMMARY 

In this paper we have outlined the results of an exten- 
sive spectroscopic campaign to map out the substructure 
of the RCS 2319+00 supercluster at z ~ 0.9. 

• From 1961 good confidence spectroscopic galaxy red- 
shifts in a 35x35 arcminute region around the three 
main component clusters of the supercluster, we have 
found 302 structure members in a redshift range of 
0.858 < z < 0.946. 

• The three clusters were found to lie at redshifts of 
0.901, 0.905, and 0.905, a separation of only 630 
kms' 1 . Combined with the < 3 Mpc transverse sepa- 
ration of each component cluster to its nearest neigh- 
bor, the three clusters of RCS 2319+00 represent the 
most compact supercluster found to date at this red- 
shift. 

• Three distinct walls in redshift space are apparent in 
the spectroscopic data, with the main supercluster 
peak surrounded by lower and higher redshift walls 
separated by ~65 Mpc (Az ~ 0.03). 

We attempt to measure accurate velocity dispersions 
from the spectroscopically confirmed galaxies for the 
three main cluster components with the aim of deter- 
mining the mass of each cluster. 
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• Cluster A: from 23 spectroscopically confirmed galax- 
ies within 1 Mpc of the cluster center, we find a ve- 
locity dispersion of 1202 ± 233 km s _1 . This results 
in a dynamical cluster mass of ~ 10 15 M©, consistent 
within errors with the lower mass estimates from the 
X-ray and weak lensing data of ~ 6 x 10 14 M Q . 

• Cluster B: there are not enough cluster members left 
within 1 Mpc of the cluster center after interloper 
clipping to calculate a reliable velocity dispersion and 
we looked to the X-ray data to corroborate our dy- 
namically derived properties. The X-ray predicted 
velocity dispersion is consistent with the measured 
spectroscopic value of 1154 ± 367 km s^ 1 and the re- 
sulting dynamical mass is consistent with the X-ray 
mass of 5.9 x 10 14 M Q . 

• Cluster C: we believe that this cluster is in a relaxed 
state and that our measured velocity dispersion of 
714 ± 180 km s^ 1 and mass estimate of 2.5 — 3.8 x 
10 14 M© is accurate. 

The multiwavelength indications of substructure in the 
supercluster field along with the close proximity of the 
three clusters in redshift and angular space, which makes 
member assignment difficult, prompted the search for 
substructure in the supercluster. 

• We performed the Dressler-Shectman test to begin 
tracing the substructure of the supercluster by look- 
ing for deviations from the global mean parameters in 
the supercluster field. We find evidence of substruc- 
ture over the broad redshift range as well as within 
each separate redshift wall, with Monte Carlo tests 
showing a very low probability that the substructure 
is false. 

• To isolate specific structure, we used a modified 
friends-of-friends algorithm in an effort to associate 
supercluster galaxies with their parent clusters or 
groups. 

• Comparing Monte Carlo simulations to the FOF re- 
covered overdensities in the real data suggests that 
groups found to have > 5 members in our real data 
are tracing legitimate structure. 

• We find evidence of possible infalling groups and ex- 
tended structure in the main clusters, with five > 5 
member groups plus the three cluster core groups, 
which we believe are tracing true structures. 

• We spectroscopically confirm the presence of a fila- 
mentary structure spanning the area between clus- 
ters A and B (group 5) previo usly detected in the 
infrared by ICoppin et al.1 ()2012h to be forming stars 
at a SFR ~ 900M Q yr _1 and building up significant 
stellar mass. The filament is detected in the DS test 
and using the FOF method. 

An investigation of the color-magnitude relation and 
color distribution of the galaxies in the different density 
environments of the supercluster was performed. 

• Of the galaxies identified as belonging to the three 
core cluster groups, 77±13% are red-sequence galax- 
ies, with blue galaxies making up only 23±7% of the 
population. 



• In the >5 member groups with the unique filamentary 
structure excluded, 43±12% of the galaxies fall on the 
red-sequence with over half of the population remain- 
ing in the blue cloud. This suggests that the group 
galaxies have undergone some processing to build up 
their red-sequence, though the remaining high frac- 
tion of blue galaxies would suggest less processing 
than the core galaxies. 

• The filamentary structure between clusters A and B 
has 8 of its 9 members below the red-sequence, with 
three distinct peaks in color: one at the bluest col- 
ors; one in the green valley; and a small peak in the 
red-sequence of the color-magnitude diagram. De- 
spite the small numbers we propose the 'green' galax- 
ies may belong t o the t ransient population identified 
by iBalogh et al.l ()2011h to be migrating to the red- 
sequence. 

We construct stacked LFs for the core and >5 member 
FOF groups from the spectroscopically weighted mem- 
bers to estimate average group richness and masses using 
two different methods: the Halo Occupation Distribution 
(HOD) method and the maxBCG richness estimators. 

• The LF derived masses from the stacked core groups 
are ~ 0.9-6x 10 14 M Q , consistent with the ~ 10 14 5 - 
10 Mq cluster masses estimates from the multi- 
wavelength data. 

• The average >5 member group masses from the two 
methods are an order of magnitude less massive than 
the core groups, at ~ 4 — 5 x 10 13 M Q . 

• Assigning cluster A as the host halo of the eventual 
merged cluster, the remaining groups account for 22% 
of the galaxies currently confirmed as supercluster 
members. Compared to studies using semi-analytical 
derived merger trees on simulated cluster data, this 
fraction is slightly higher than the fraction seen at this 
epoch in massive, present-day clusters but not unreal- 
istic given that all of the studies found a large fraction 
of ~ 30% of surviving galaxies in present day clus- 
ters spent some time as members of ~ 10 13 h^ 1 M© 
groups. 

The mapped structure of the RCS 2319+00 superclus- 
ter defined in this paper provides an abundance of differ- 
ent density environments in which to study the effect of 
the environment on galaxy evolution, from the isolated 
supercluster field galaxies, along infalling groups and fil- 
aments to the dense cluster cores. Our vast catalog of 
multiwavelength data on this field will be used in con- 
junction with the spectroscopic catalog defined in this 
paper to study the environmental effects on our mem- 
ber galaxies (e.g. a study of the star-formation and stel- 
lar mass in RCS 2319+00; Faloon et al. 2012, in prepara- 
tion). Though again we caution that the structure de- 
fined here is a lower limit on the possible groups present 
in the RCS 2319+00 supercluster field and 'isolated' su- 
percluster galaxies may not actually be isolated given 
further spectroscopy. The addition of high-quality, mul- 
tiwavelength photometric redshifts and supplementary 
spectroscopy will aid in limiting the mis-identification 
of the different density regions within the supercluster 
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and add to our understanding of this unique, compact, 
high-redshift supercluster. 
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